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Porous media, Fast diffusion equations and the existence of 
global weak solution for the quasi-solution of compressible 

Navier-Stokes equations 

Boris Haspot * 

^-^ . Abstract 

CN . In [T71 [T51 [201 im 122] , we have developed a new tool called quasi solutions which 

^H ' approximate in some sense the compressible Navier-Stokes equation. In particular it 

M-i, allows us to obtain global strong solution for the compressible Navier-Stokes equa- 

■^>( ' tions with large initial data on the irrotational part of the velocity {large in the sense 

^^ , that the smallness assumption is subcritical in terms of scaling, it turns out that in 

this framework we are able to obtain a family of large initial data in the energy space 
in dimension N = 2). In this paper we are interested in proving the result anounced 

^ ' in [53] concerning the existence of global weak solution for the quasi-solutions, we also 

observe that for some choice of initial data (irrotationnal) the quasi solutions verify 
the porous media, the heat equation or the fast diffusion equations in function of the 

r^ , structure of the viscosity coefficients. In particular it implies that it exists classical 

C^ ' quasi-solutions in the sense that they are C°° on (0, T) x M^ for any T > 0. Finally 

we show the convergence of the global weak solution of compressible Navier-Stokes 
equations to the quasi solutions in the case of a vanishing pressure limit process. In 
particular we show that for highly compressible equations the speed of propagation 
of the density is quasi finite when the viscosity corresponds to /i(p) = p" with a > 1. 

fvl . Furthermore the density is not far from converging asymptotically in time to the 

f^ ' Barrenblatt solution of mass the initial density pQ. 

in 

t:;^' ■ 1 Introduction 

o' 

en I The motion of a general barotropic compressible fluid is described by the following system: 



> 



X 



' dtp + div{pu) = 0, 

< dtipu) + div{pu ^u)- div(2/i(p)L>(n)) - V{X{p)dwu) + VP{p) = pf, (1.1) 

^ {p,u)/t=o = {po,uo)- 

Here u = u{t, x) G M stands for the velocity field, p = p(t, x) G M"*" is the density 
and D{u) = ^{^u -|-* Vu) the strain tensor. The pressure P is such that P{p) = ap^ 
with 7 > 1 and a > 0. We denote by p.[p) and \[p) the two-Lame viscosity coefficients 
depending on the density and satisfying: 

p{p) > 2p{p) + NX{p) > 0. (1.2) 
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Throughout the paper, we assume that the space variable x E M . In this article, we 
are going to investigate the existence of global weak quasi solutions for the system (jl.l|) , 
a notion which has been introduced in |21l [22l [T71 [T8] in order to prove the existence of 
global strong solution with large initial data on the rotational and irrotational part of 
the velocity for the scaling of the equation (we refer to the remark H] for the definition 
of this notion of scaling when we assume that the initial density is far away from the 
vacuum) . 

Before entering in the heart of the topic and defining the notion of quasi-solutions, we 
would like to explain more in details the notion of invariance by scaling for compressible 
Navier-Stokes equations when we suppose that the density po is in a Besov space. Up our 
knowledge these results are quite new and shall allow to understand the relation between 
the quasi solutions and the so called Barrenblatt solution for the porous media equation 
in terms of invariance by scaling. 

Scaling of the equations 

A natural way to understand the equations of fluid mechanics correspond to search self 
similar solution, it means that there is a scaling of the variables after which the system 
become stationary solutions. Precisely it holds when we set: 

p{t,x) = r^F{xt-^), 

u{t,x) = t-''^G{xt-^) 

The exponents a, ai and /3 are called similarity exponents, and functions F and G are 
the self similar profiles. In particular a and ai are the density contraction rate and /3 
the space expansion rate. In the sequel we shall assume that /i(p) = fip^ and A(p) = Xp^ 
with 6 > 0. Simple calculus give when we set r] = xt~^: 

dtp{t, x) = -t-'^-^aFir]) + /3VF(77) • 7?), 
div(H = t^"^"i-'^div(GF)(r/). 



Next we have: 



and finally: 



pdtu = -t-"-"i-i(aiF(r/)G(??) + /3r/ • VG(r?)F(7?)), 
pu-Vu = t-''-^'''-'^F{ri)G{r]) ■ VG{rj), 
2ndw{p^Du) = 2/it-^"-"i-2/5div(F''(??)Z?G(r/)), 
V(/divn) = t-^^'^^'^'^V{F\7])diYG{r])). 



In order to ensure the existence of solution for the svstem (|l.ip under the specific form 
introduced in ()1.3p we need to assume that: 

= a + 1, 

: a + 2ai + /3, 
■-6a + ai + 2/3, 
= 07 + /3, 




which is equivalent to the following system: 



ai + /3 = 1, 
((9-l)a + 2/3 = 1, 
a(7-l) + /3 = ai + l. 

The solution of the previous system is: 

-1 1-7 „ 26* -7-1 , , 

With this choice on the parameter a, ai and /3 we finaly get the profile equation: 

aF{r]) + l3VF{r]) ■ r] - dw{GF){r]) = 0, 

aiF{rj)G{7]) + /3r/ • VG(r/)F(7?) - F{r])G{r]) ■ VG{r]) + 2/idiv(F^(7?)Z)G(r/)) 

+V(F^(r/)divG(ry)) - V(aF(r/)T) = 0. 



Remark 1 A first remark consists in observing that there is no scaling invariance when 

An other way to express this scaling invariance corresponds to consider a classical solution 
{p, u) of the system (|1.1|) and to check that the family: 

{pi,ui){t,x) = {rp{it,i^x),np{it,i^x), 

is a solution of (jl.ip for any / G M when a, ai and /3 verify ()1.4p . 

Remark 2 We shall say that a functional space E embedded in S (R^) x [S (M^)))^ 
is a critical space for il.l\) if the associated norm is invariant under the transformation 
{po,Uo) — )• {{po)i, iuo)i) for any / G M. In particular we observe that: 

JV I 2 JV 1 I 2_ s-1 



D 2 ^ 29-7-1 ^ r2 -^"^^29-7-1 
-D2,oo ^ -02,00 ' 

verifies a such property ( we refer to J^ for the definition of Besov space) . 

Remark 3 Let us point out that when p[p) = p,p and A(p) = with 7 > 1 (the case 

of the shallow-water equation) then -62^00 ^~^ ^ -^2^00 ^-^ ^ critical space invariant by the 
scaling of the equation. In particular let us mention that it is the same scaling invariance 
for the initial velocity uq than for the incompressible Navier-Stokes equations. 

JV__2_ JV -. I 2 

In the case of constant viscosity coefficients we observe that i?2^oo ^^^ ^ -^2^00 ^^^ ^■^ ^ 
critical space for the scaling of the system il.l\) . 

Remark 4 To finish we would like to mention as it has been observed in f^ that there is 
no invariance by scaling when we wish to work with density far from the vacuum,, typically 
po = 1 + qq with qq in a Besov space. However in 19] Danchin makes abstraction of the 



pressure term and define a other notion of critical space for the compressible Navier- 
Stokes equation when the initial density is far away from the vacuum. More precisely we 
can remark that lil.l\) is invariant by the transformation: 

{poix),uo{x)) -> {pq{Ix),Iuq{Ix)) {p{t, x),u{t, x)) -^ {p{l'^t, Ix), lu{l^t, Ix)) 

up to a change of the pressure law P into PP. In particular B21 x B2 i is norm 
invariant by the previous transformation and is critical for the initial data {qo,uo). This 
notion of critical space seems well adapted to the case of initial density far away from the 

vacuum and is also relevant at least for the initial velocity since B21 is also critical for 
the incompressible Navier-Stokes equations. Let us recall that with this notion of critical 
space we have proved the existence of global strong solution with large initial data on the 
irrotational and rotational part of the velocity (see \11\ \1^ \20f ) by involving the notion 
of quasi- solutions verifying regularizing effects. 

Before coming back on this notion of scaling in order to make some link between self 
similar solutions for the porous media equation and the behavior of the quasi-solutions, 
we would like now to give precise assumptions on the viscosity coefficients with which we 
are going to work. 

Condition on fi{p) and X{p) 

In this paper we are interested in extending the notion of quasi-solution developed in 
|2 H I22 1 [T7 1 ll8] (where only shallow water coefficients were considered) for general viscosity 
coefficients following the algebraic equality discovered by Bresch and Desjardin in |6ii7j: 

X{p) = 2pfi'{p)-2fi{p). (1.5) 

In the sequel we shall deal with the function (p{p) defining by (p (p) = ^ ^'^' and the 

function f{p) by / (p) = y/p^p (p). Let us mention that the equality ()1.5p implies that the 
viscosity coefficients are degenerated inasmuch as it imposes that /i(0) = A(0) = 0. With 
this choice of viscosity coefficients Bresch and Desjardin have obtained a remarkable new 
entropy for compressible Navier-Stokes equations (jl.ip providing a Ly'(L^(M^)) control 
for any T > on the gradient of the density (more precisely on yfpV'^{p)). In particular it 
allows them to prove the existence of global weak solution for a specific choice of pressure, 
more precisely what they describe as a cold pressure. Compared with the case of viscosity 
coefficient the pressure term is quite simple to deal with by using Sobolev embedding 
since we have uniform estimate on yjp\/ip{p) in L^(L^(M^); however a new difficulty 
appears coming from the degenerescence of the viscosity coefficient. Indeed we lose the 
control of Vu G -^^((0, T) x M^) what makes delicate the compactness study of the term 
pu®u (in particular via the classical energy estimates we have only a convergence in the 
sense of the measure) due to the existence of vacuum. In order to overcome this difficulty 
Mellet and Vasseur in [29j obtained new entropy on the velocity which furnish them a 
gain of integrability on the velocity. With this new ingredient they are able in [29j to 
prove the stability of the global weak solution for compressible Navier Stokes equations 
with such viscosity coefficients and with classical 7 law pressure P{p) = ap^ with a > 



and 7 > 1. 

We are going to detail here the assumptions on the viscosity coefficients which ahow to 
Mehet and Vasseur to obtain additional informations on the integrability of the velocity 
and we are going even to relax their hypothesis (it will be crucial in order to consider 
in the sequel quasi solution verifying porous media equations). We shall suppose the 
following inequalities on /x and A, let ui £ (0, 1) and f2 > such that: 

|A'(p)| < -fi'ip), 

vi (1.6) 

Vl^x{p) < 2p{p) + N\{p) < U2p{p)- 

Remark 5 // we assume that p[p) = pp" with a > then the relation il.5\) gives: 

X{p) = 2{a-l)pp", (1.7) 

and: 

2p{p) + N\{p) = 2(1 + N{a - l))p{p). (1.8) 

In this situation we have vi = U2 = 2(1 + N{a — 1)). 

Following [29j let us briefly make some comments on the conditions ()1.6p . 

Remark 6 The condition < |j.6]) is crucial in order to obtain the estimates \3.S9\I and 
li3.38\) . In particular we can observe that the second condition in \1.6\) is similar to the 
classical assumption on the Lame coefficient 2p{p) + NX{p) > when p{p) = pp" and 
\{p) verifies ( fi.5|) . 

Remark 7 The lower estimate in the second inequality in M.(?^) is trivial when \[p) > 0, 
while the upper estimate is trivial when X{p) < 0. Together this provides: 

\\ip)\<Cp{p) Vp>0. 

This inequality and the first inequality of \1.6\) will be crucial for estimating the limit of 
V {X{pn)divun) . 

Remark 8 Condition lil.6\) and m.5\) implies that: 

N-i + iq- p'{p) iv-i + if 

^ < ^-f^ < ^ Vp > 0. 

Np - p{p) - Np ^ 

It yields: 

j Cp^-^+5^ < p{p) < Cp^-^-^^ Vp > 1, 

\ Cp^-^+5^ < p{p) < C7pi-^+5^ Vp < 1, 

We can now recall briefly the definition of the quasi solutions introduced in \n\ [T8\ \ZT\ [22] 
which doughly spaaking are solutions of the compressible Navier-Stokes equations (jl.ip 
where we have cancelled out the pressure term (in the sequel we shall give a more accurate 
definition) . 



Definition 1.1 We say that {p,u) is a quasi solution if {p,u) verifies in distribution 

sense: 

( d 

—p + div{pu) = 0, 

— (pn) + dw{pu (S)u)- div(2/i(p) Du) - V(A(p)divu) = 0, ^^'^^^ 

_ (yO,n) t=o = {po,uo) 

As we explained previously this notion of quasi solution is interesting inasmuch as it 
allows to exhibit large initial velocity on the irrotational part in the sense of the scaling 
of the remark [J] (in particular we assume no vacuum on the initial density) providing 
global strong solution for compressible Navier-Stokes equation (see [T71 [181 EO] for more 
details). This result is based on a strange phenomena on the quasi solution since these 
last one verify regularizing effects allowing to neglect in terms of scaling the pressure 
term in high frequencies. An other way to express the things is that the quasi solutions 
preserves a structure of irrotationality of the system when we choose irrotational initial 
data (it will be the case in the sequel), we shall say that the quasi solution, typically 
u = —'V<p{p) is purely compressible. For other results on the existence of strong solution 
with critical initial data for variable viscosity coefficients we refer to |8l [T5l [T6] . 

We now are going to investigate the existence of such quasi solution for the viscosity 
coefficients verifying p.5|) when the initial data is assumed to be close from the vacuum, 
typically po ^ L^(M^). More precisely as in [T71 [181 [2D] we are going to search in a first 
time irrotational solution under the form u{t,x) = 'Vc{t,x) for the system (|1.1U|) . Let us 
assume now to simplify that: 

l^(p) = pp"" with a > and \{p) = 2(a - l)pp°', (1.11) 

with a > 1 — jy in order to ensure the relation 2p{p) + N\{p) > 0. We observe here 
that p{p) and X{p) verify the relation (jl.Sp . In this case we will verify (see the theorem 
II. ip that at least for suitable initial data on po then it exists a explicit solution to the 
problem (jl.lOp written under the form (p, — V(/?(p)) with p verifying the porous media or 
the fast diffusion equation when a / 1: 

- °' (1.12) 

In a very surprising way it means that the quasi solutions are directly related to the 
porous or the fast diffusion equations. Moreover we will show that when we work with 
highly compressible Navier-Stokes equations (which correspond to the case where a goes 
to with -P(p) = ap'^), then the properties of porous media or fast diffusion equations 
are more or less preserved for the solution of (jl.ip when a is small. Before giving more 
details on your results, let us give few words on the porous and fast diffusion equations 
for the reader which are not so familiar with these equations. 

Porous media and fast diffusion equations 

Let us consider the equation (J1.12p when 2/i = 1 to simplify the notations; the case 
a > 1 (the porous media equations) arises as a model of slow diffusion of a gas inside a 




porous container. Unlike the heat equation a = 1, this equation exhibits finite speed of 
propagation in the sense that solutions associated to compactly supported initial data 
remain compactly supported in space variable at all times (see )33j and [1]). When 
< a < 1, the opposite happens. Infinite speed of propagation occurs and solutions 
may even vanish in finite time. This problem is usually referred to as the fast diffusion 
equation. 

Let us recall that there exists a theory of global unique solution for initial data po ^ 
L^(]R^) (see the section [2] for some reminders). M. Pierre in [31] has extend this last one 
in obtaining the existence of unique global weak solution with bounded Radon measure 
as initial data. Let us mention also that the porous media equations are invariant by 
scaling, more precisely we can introduce a notion of self similarity (for more informations 
we refer to the chapter 16 of |33]). The notion of scaling consists in searching some 
solutions under the following form p{t, x) = t~'^F{^), with 7 and f3 to be determined. In 
our case 7 and (3 have the form: 7(0 — 1) + 2/3 = 1, and F verifies the following equation: 

AF° + I37]VF + 7F = 0. 

In this case we said that p is a self similar solution of type I or a forward self similar 
solution. In particular it exists self similar solution such that the initial data is a Dirac 
mass (as in the theorem of M. Pierre in j31j) and such that for t > this solution 
conserves a constant mass. This is the so-called Barrenblatt solutions (here a > 1) that 
we can write under the following form: 

U^{t,x) = t-^-F{^) with F{x) = {C-^-^^^\xf)t~' 

with C > and 71 = ^y/ _-^^ , 2 , P = 7vTa^TTP2' Here we have the conservation of the 
mass / Umit, x)dx = m with m depending on C and the initial data corresponds to the 
Dirac mass m6o. Similarly when rric < a < 1 with rric = max(0, — ^) it exists also 
Barrenblatt solutions defined as follows: 

Um{t,x) = t-^^F{xt-'^) with F{x) = {C + Ki\x\^)^ , 

with Ki = 2Na • ^^ recall that asymptotically in time all the global weak solution with 
L^ initial data converges to a Barrenblatt solution determined by his mass ||tio||Li (we 
refer to Friedman and Kamin [14j, Vazquez and Kamin [241 [25] and Dolbeaut and Del 
Pino [12] ) . As we mentioned previously in the case of fast diffusion equation < a < 1 , 
infinite propagation occurs and solution may even vanish in finite time when a is in the 
interval (0, rric) with rric = max(0, -^—)- In particular it implies a lost of the initial mass 
when pq is in L^ (it implies also a lost of the regularity of the solution). We refer to 
[34 ] theorem 5.7 for a necessary condition of extinction, in particular the initial data 



belongs in an appropriate Marcinkewitz space Mp*(M ). Let us finished this subsection 
by mentioning that we shall recall more results on the porous and the fast diffusion 
equations in the section [21 

1.1 Main results 

Let us now give your first result describing the link between quasi-solutions and the 
solutions of ()1.12p . Finally we obtain the following theorems. 



Theorem 1.1 Let ^{p) = fip'^ with a > 1 — j^ and X{p) verifying the relation il.5\) . 
Let pq G L^(R^) with pQ > and continuous and uq = —'S/ip{pQ). Then it exists a global 
weak solution solution of the system \1.10\) of the form {p,u = —Viplp)) with {p,u) 
belonging in C°°((0,+oo) x M^) n C([0,+oo] x M^) and solving the following system 
almost everywhere : 

{ dtp - 2Ap{p) = 0, 



(1.13) 



Furthermore we have: 



lim \\p{t)-U^{t)hi(^^N^ = 0. (1.14) 



Convergence holds also in L°° norm: 

lim t^\\p{t)-U„,{t)\\Lo.^^N)=0, (1.15) 



t— >+c« 



with /3 = yy/ _i-,_i_2 ^^'^ ^m the Barrenblatt of mass m = ||po||li(E^)- -^'^'^ every p G 



N{a-l)+2 

(1, +oo) we have the following regularizing effect, p{t, •) belongs in LP{M.^) and: 



with Gp - (Ar(„„i)+2)p "'^d ttp - ^]sf(a-l)+2)p ■ 

Remark 9 Let us point out that we could have also global strong solution for more 
general viscosity coefficients with p. verifying the same conditions than the subsection 
\2.1.1\) . We refer to the section\^for more details in this situation. 

Remark 10 We shall remark in the proof of this theorem that any solution of p.l2p 
such that p is in C^((0, +oo) x R ) is a solution of /IL10\) almost everywhere. 

Remark 11 We can observe as in jildf that if we consider the compressible Navier-Stokes 
equation with a friction term apu and a pressure term of the form 2pap°' then the solution 
of the previous theorem \1.1\ verify also a such system. 

Remark 12 Let us mention that when a is in the interval (0, rric) with nic = inax(0, -^—) 
(then this contradicts the Lame condition on the viscosity coefficients ) then it can ap- 
pears a phenomena of extinction of the solution in finite time, in particular it implies a 
loss of the initial mass when po is in L^. A typical example is the solution: 



;^)i^ and u{t,x) = -^V(Ji: 
xr a — 1 1 — t 



with Cq " = 2[N — jr^)- In particular we observe a blow-up behavior of u at time T. 

We are going to give a general definition of global weak solution for the quasi solutions in 
the spirit of [2U] including the case where the initial velocity is not necessary irrotational. 



Definition 1.2 We say that {p,u) is a global weak quasi solution if {p,u) verifies in 
distribution sense: 

( d 

—p + div{pu) =0, 

— (pn) + div(pn (g) u) - dw{2p{p) Du) - V(A(p)divu) = 0, ^^'^^^ 

_ (p,n) t=o = iPo,uo). 

More precisely {p, u) is a weak solution of il.lU\) on [0, T] x M^ 

P/t=Q = Po > 0, pn/j=o = "T-o- (1-17) 

with: 



po e L^M^), ^V^ipo) e l2(m^), po > 0, 



(1.18) 



/9 G L5?(Li(M^), ^pV99(/9) G L~(L^(R^)), ^u G L~(L^(R^)), 



• VKp) Vw G i'((0,r) X M^), ^|n| VRTTH^ G Lf{L\R^)). 

with p>0 and {p,^u) satisfying in distribution sense on [0,T] x R^; 

f dtp + div(^^u) = 0, 
\p(0,x) = po{x). 

and if the following equality holds for allip{t,x) smooth test function with compact support 
such that (p{T, •) = 0; 

{pu)o ■ ¥?(0, ■)dx + / / ^{y/pu)dtip + ^/pu ® ^u : Vipdx 

Jo Jrn (i-iyj 

- < 2p{p) Du, \7(f > - < X{p) divu, dwip >= 0, 
where we give sense to the diffusion terms by rewriting him according to ^/p and y/pu: 

— —{y/puj)dii(pjdxdt— I 2{y/puj)p {p)diy/pdiipjdxdt 

-{y/puj)dji(pjdxdt — / 2{^Ui)iJi {p)dj^diipjdxdt 



/\ ( \ r 

— —{,/pui)djHPjdxdt— / 2{^/pui)X {p)diy/pdjipjdxdt 

We assume also that p and A verify the conditions il.5\} and il.6\) . 

We obtain now a general result concerning the stability of the global weak solution for 
system (|1.10p and a result of existence of global weak solution for general initial data of 
the form (/jq, —'V(p{pq)) (in particularly po is not assumed only strictly positive). 



Theorem 1.2 Assume that i^{p) and A(p) are two regular function of p verifying lll.5\) 
and lll.6\) . Furthermore we shall set g{x) = ^^^ and we assume that g is bijective and 

that g~^ is continuous on (0, +cxd). When 2 + N < ui, we assume in addition that g and 
g are increasing on (0, +cxd). 

Let {pn,Un) be a sequence of global weak solutions of system \1.24^ satisfying entropy 
inequalities i3.36\) . ^3.37^ and /i3. 38]) with initial data: 

{Pn)/t=0 = Poi^) <^nd {pnUn)/t=0 = Po'^oi^) 

with Pq and Uq such that: 

Po > 0, p^ ^ po tn Li(E^), p^u'S -^ pouo m L^R^"), (1.20) 

and satisfying the following bounds (with C constant independent on n): 

[ pS^<C, / v^|Vv.(pS)pdx < C (1.21) 



and: 



I 

Jr' 



1 _L |7/"|2 

pg ^' °' ln(l + |u[f |2)dx < C. (1.22) 



Then, up to a subsequence, (pn, y/fhiun) converges strongly to a global weak solution 
{p, -./pu) of {1.24^ satisfying entropy inequalities 113.36]) . 1^3.37] ) and i3.38]) . 
Furthermore the density pn converges strongly to p in C([0, T], Lj^"(]R^)) with < a < 
z/i when N = 3 and in C{[0,T],Ll^^{W^)) for any q>l when N = 2; ^/p^Un converges 
strongly in L^{0,T,Lf^^) to ^fpu and the momentum nin = PnUn converges strongly in 
L^{0,T,Ll^{R^)), for any T > 0. 

If we assume moreover that {pq,uq) verify the initial condition of the definition M.SX that 

uq = —Vip{p) with p{p) = pp" and that y^|uo| ^ belongs in L^(R^) forp large enough, 
then it exists a global weak solution (p, u) of the system \1.24^ where p is also the unique 
solution of the system il.l2]) (see the theorem 2.4] for the existence of a unique global 
strong solution for the equation \1.1^) with a L initial data). 

Remark 13 Let us mention that the additional technical assumption on the viscosity 
p remains quite natural since they are verified in the standard case p{p) = pp" with 
a > 1 — jy . In particular we remark that this result extend the analysis of \29^ to general 
viscosity coefficient, in particular we do not suppose that p (p) > c > 0. 

Remark 14 Let us emphasize that the condition M.(?(l implies that we exclude the case 
of fast diffusion equation with ^ < a < 1 — jj, in particular it forbids any phenomen of 
extinction and loss of mass ||p(i, •)||li(K^)- 

Remark 15 Let us mention that our second result of existence of global weak solution 
when uo = — V(/?(po) can be applied to the Barrenblatt solution by choosing po = UmiT, •) 
with r > 0. In this case we have then a free boundary problem with (p, u) C°° when p > 
and {p,u) = (0,0) when p = 0. In particular when p(p) = pp" with a > 1 we observe 
that the propagation of the support of the Barrenblatt is finite. 

10 



More generally when we choose a initial density with compact support the support of 
the density remains hounded along the time when we deal with the case fi{p) = /ip" with 
a > 1. It means that quasi- solutions in this case have the same properties than the porous 
media solutions. We are related to a free boundary problem with on a side the solution 
which is C°° and on the other side the solution which is identically null. In particular it 
implies that we can not hope uniqueness as the velocity can take any value on the vacuum 
set, that is why it is natural to consider the momentum unknown ^Jpu as it is the case 
in the previous theorem. 

Remark 16 Let us emphasize on the fact that the problem of the existence of global weak 
solutions remains open in the general case (it means when uq is different from —'V(p{pQ)). 
Indeed in the previous theorem we prove the stability of the global weak solutions, however 
it seems quite complicated to construct approximate global weak solution which verify 
uniformly in n all the entropies. We have the same problem in the case of the compressible 
Navier Stokes problem where the problem remains also open. 

We are going to prove now the convergence of the global weak solution of the compressible 
Navier-Stokes equations to the quasi solutions when we consider a vanishing pressure 
process. More precisely let us consider the highly compressible Navier-Stokes system 
with e going to 0: 

' dtPe + div(/>eUe) = 0, 

dt{p,u,) + dw{p,u, ® u,) - div{2p{p,)D{u,)) - V{X{p,)dwu,) + eVp] = 0, (1.23) 

, iPe,Ue)/t=0 = {PQ,Uo)- 

In the literature we find many result on the incompressible limit which corresponds to 
take e = ^z with r/ the Mach number going to 0. For such results in the framework of 
the global weak solutions for the ill-prepared data we refer to the pioneering papers of 
Desjardins, Grenier, Lions and Masmoudi [lOl [TTl [271 [28] . Indeed in this last situation we 
can observe at least heuristically that the density p^ converges to a constant 1 when we 
are working with initial density of the form po,^ = 1 + ^5o,e with go,e uniformly bounded 
in appropriate space. By considering the mass equation, roughly speaking we can deduce 
that the limit solution of u^ is incompressible. The main difficulty compared with the 
well prepared case corresponds to deal with the acoustic waves, in particular in order to 
overcome such difficulty the authors use Strichartz estimates. 

In our case we are going to deal with the opposite situation when the solutions are highly 
compressible and converge to quasi-solution which are in some sense purely compressible 
since irrotationnal. In the case of constant viscosity it appears impossible to pass in the 
limit when e goes to since we lose the L^(L'^(M^)) estimate on p^ coming from the 
pressure term (we conserve only the L^ conservation of the mass which is not sufficient 
to pass to the limit since it does not provide enough compactness information). In our 
case with viscosity coefficients verifying the relation (|1.5|) we know that we have uniform 
estimate on y/p^Vip{p^) in L^(L^(R^)) providing of the entropy ()3.42p (see [6]) which 
will be sufficient in term of compactness to pass to the limit when e goes to 0. Before 
giving our main result in this spirit, let us give a definition of global weak solution for 
the system ()1.23p in the spirit of [29] : 
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Definition 1.3 We say that {p,u) is a global weak solution of lil.l\) if {p,u) verifies in 
distribution sense: 

( d 

—p + dw{pu) = 0, 

— ipu) + div(pn n) - div(2^(/)) Du) - V(A(/))divu) + VF{p) = 0, ^^-^^^ 

ip,u) t=o = ipo,uo), 

with P{p) = ap"^ , 7 > 1. More precisely {p,u) is a weak solution of Hl-l]) on [0,T] x R 

P/t=o = Po > 0, pn/j=o = "T-o- (1-25) 

with: 



po G L'{R'') n L-'iR''), ^V^ipo) G L'{R''), po > 0, 
if 



,|(i + 7mTTkF))gl'(m^) 



(1.26) 



p G L^{L^R^), ,/pVipip) G L^(L2(M^)), ^n G L^{L^{R^)) 
y^VuGL2((o,r) xM^), ^|n|v^ln(l + |u|2) e LfiL"^ 



with p> and (p, ^pu) satisfying in distribution sense on [0,T] x R^ : 

J dtp + div{^y/pu) = 0, 
\piO,x) = po{x). 

and if the following equality holds for allLp{t,x) smooth test function with compact support 
such that ip(T, •) = 0; 

{pu)o ■ ip{0, ■)dx + / / ^{y/pu)dt<p + y/pu ® y/pu : Vifdx 

Jo Jm (-^27) 



rT 



+ a / p'^diY(p{s, x)dx— < 2p{p) Du, Vy? > — < \{p) divn, div(^ >= 0, 

where we give sense to the diffusion terms as in the definition MJA We assume also that 
p and A verify the conditions M.5\) and \1.0\) . 

Theorem 1.3 Let^ > land {pQ,uo) verifies the conditions of definition \1.3[ Let us take 
the following assumptions on j, vi and 1^2 •' 

• iV = 3 
• vi> 2: 
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• < i/i < 2; 



N 



6 12 ' 2 - z^i ■' 

6 12 ' 6 12 -^ 



1 + |<7 ^/iV = 2. 



T/ien under these conditions we have as in the theorem \l.S\ the stability of the global weak 
solutions for the system il.l]) . 

AssuTne that there exists global weak solution {p^,u^) verifying the definition of 129] with 
the conditions M.5\) and M.0\l on fi{p) and X{p). Then {p^,Ue) converges in distribution 
sense to a global weak quasi- solution {p, u) of the system {1-24^ in the sense of the defini- 
tion \1.2[ Furthermore the density p^ converges strongly to p in C([0,T],L^^'^(]R^)) with 
< a < ui when N = 3 and in C([0, T], L^^^(M )) for any q > I when N = 2; ^/p^u^ con- 
verges strongly in L^(0, T, L^^^) to y^ti and the momentum m^ = p^^u^^ converges strongly 
in Li(0,T,Li^^(M^)), for any T > 

Remark 17 Let us mention that the technical assumption on ui, U2 andj are important 
in order to ensure a uniform gain of integrability on the velocity u^ (as in the paper of 
Mellet and Vasseur f29^). more precisely we will see that we have a competition between 
the pressure and the viscosity. 

Remark 18 We shall also emphasize on an important question which remains open; in- 
deed when p{p) = p.p"' and uq = —Vip{pQ) the limit solution {p,u) of {pe,Ue) when e goes 
to is a quasi solution of il.l3\) . However it is not clear how to prove the uniqueness 
of the quasi solutions in the class of the solutions giving by the definition \1.2\ and in 
particular to show that this quasi solution is solution of the porous media equation for the 
density p when a > 1. In the following corollary, we shall give properties of the solutions 
of M.2'J\) when we assume the uniqueness of the quasi solutions. 

Let us mention that in some case if we assume that the solution of the porous media 
is enough regular (typically the case of some solutions with initial data po = Um{T,-) 
with T > and Um. o- suitable Barrenblatt solution) the uniqueness consists in prov- 
ing a weak-strong uniqueness theorem since we can assume that the solution (Um{T + 
r, •), -V(/?(f7m(T + T, •))) is strong. 

Remark 19 The second important remark consists to point out the fact that the question 
of global weak solution for the system \1.2'J\l remains open, indeed Mellet and Vasseur 
have proved the stability of the global weak solution in 129^ . However it seems not so 
easy to construct a regular sequel {pn,Un) approximating \1.2'J\) (typically by a Friedrich 
process) and verifying uniformly all the entropies of f2^ i3.41\ ), i3.42\ ) and 1^3^43^. 

Corollary 1 Let 7 > 1, i^i and V2 with the hypothesis of theorem li.3l Let pQ and uq 
verifying the hypothesis of the theorem \L3\ with uq = — V</j(po) o-iT'd p-ip) = p-p°' , A(p) 
verifying also the hypothesis of theorem IJ.31 We assume here that there exists a unique 
quasi solution of system \l-24^ with such initial data, in particular the density of this 
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quasi solution verifies hl.l2\l . Then as in theorem \1.3\ (p.. uA converges in distribution 
sense to a global weak quasi- solution {p, u) of the system ^1.24\ ) such that p is solution of 
\1.12\) . Moreover: 



If a > 1, po ^ L°°(M ) and the support of po is compact then we have for T > p^ 
converges strongly to p in C([0, T],L^{R^) D LP{R^)) with l<p<l + uiifN = 3 
and in C([0,r],LP(R^)) for any p > 1. 
For all rj > it exists eo > such that for all < e < eo \/t £ [0, T] we have: 

p,{t,-)=p{t,-)+fe{t,-), 

with p{t, •) with compact support for any t G [0, T] (it exists C > such that 
suppp{t, •) < Ct^ with 13 = jjj^^rT)+2) "^'^ ^^ ^"^^ ll/^(*' OIIlicr^) < ?? Vt G [0,r]. 

In particular for all rj > it exists eo > such that for all < e < eq we have: 

p,{t,-)=p{t,-)+fe{t,-), 

with p{t, •) verifying for p as above: 

'^^^^ ^p = (lila-ilXl^p ""'^ "p = (jv(^-i~)+2)p ' "-^^ '^^ ^"^e ||/,(t, OIIlpcr^^) < n for 
allt G [0,T]. 

If oi> \ — -^, for all rj > for all compact K it exists T > 0, it exists eo > such 
that for < e < eo we have: 

\\Peit,-)-U^{t,-)\\LHK)<r] VtG [T,2T]. 



Remark 20 It is very surprising to observe that for e small enough p^ is subjected to a 
type of damping effect in L^(M ) with p choose as above up to a small remainder term 
in LP(M ). Let us point out that this effect is similar to the dispersion property for the 
Schrodingier or the waves equations. In fU^ \15f we observe a damping effect on the density 
due to the role of the pressure, but in our case the pressure tends to disappear. As for 
the porous media equation this effect seems purely non linear and is exhibited because the 
particular structure of the viscosity coefficients. 

Remark 21 Under the hypothesis of uniqueness of the quasi solution when uq = —'V(p{pQ) 
we show that the solution of the highly compressible Navier-Stokes equation are not so far 
to have a finite speed of propagation when we take a initial density with compact support. 
Indeed this is the case modulo a perturbation f^ of small L^ norm. Similarly modulo 
this hypothesis of uniqueness we expect a asymptotic convergence of p^ to the Barrenblatt 
solution of L^ norm H/OollLimiV) modulo a small perturbation. 

The paper is structured in the following way: in section [2] we recall some important 
results on the porous media and the fast diffusions equations. In section [3] we adapted 
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the entropy of [6j and [29j to the case of the quasi-solutions. In section SI we give a 
few notation, some compactness results and briefly introduce the basic Fourier analysis 
techniques needed to prove our result. In section S] we prove theorem II .11 and in section 
[5] we show the theorem 11.21 In section 15.11 we conclude with the proof of the theorem 
11.31 and the forollary [TJ An appendix is postponed in order to prove rigorously some 
technical lemmas. 

2 Important results on the porous media equation 

for the sake of completeness for the reader which are not familiar with the porous and 
the fast diffusion equations, we are going to recall in this section some essential results 
on the porous media and the fast diffusion equations. The majority of them are directly 
issue from the excellent book [33], [33| from Vazquez. In this part in order to simplify 
the problem we shall only consider the following equation with a > 1 — ;^: 

\ p(0, •) = /oo > 0. 
Let us start with the case where a > 1 

2.1 Porous Media, a > 1 

In the sequel we shall set Q = (0,+oo) x R^. Let us recall the notion of global strong 
solution for the equation ()1.12p of the porous medium equation (a > 1) (see [33] chapter 
9 for more details). 

Definition 2.4 We say that a function p G C([0, +oo), L^(M^)) positive is a strong L^ 
solution of problem i2. 28\) if: 

• /9" e Ll^{0,+oo,L^{R^)) andpt,Ap^ G L[^^((0,+oo) x 

• pt = /uAp" in distribution sense. 

• u{t) -^ pq as t ^ in L^ 



Let us mention (see [33j pi 97) that we have the following theorem of global strong 
solution, we are going to give a sketch of the proof of the existence and uniqueness of 
the L^ solution (which is perfectly detailed in |33j Chapter 6 and 9). Indeed it will be 
important to understand this point for the proof of the last part of theorem 11.21 

Theorem 2.4 Let a > 1 For every non-negative function po G L^(M^) n L°^(IR^) there 

exists a unique global strong solution p > of l[2.28]} . Moreover, dtp £ ^fociQ) f'^''" 

1 < P < ^ and: 

dtP>- , \.^ in V\Q), 
[a — \)t 

ii« u Ml ^ 2||/)o||l1(k'^) 
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Let pi and p2 be two strong solutions of 112. 28\) in (0, T) x M then for every < r < t 
\\{piit,-) -P2(i, •)) + ||Li(RiV) < ll(pi(^, •)-P2(r, •)) + ||li(k^')- (2-29) 

// pi and p2 are two strong solution with initial data poi o-i^d po2 with poi < Po2 in M.^ , 
then pi < p2 almost everywhere in (0, +cxd) x M . 

Proof: The proof is decomposed in three steps. Let us recah that in the sequel the 
initial density is always positive. 

Energy global weak solution on a bounded domain Q 

Let us mention that a global weak solution in our setting is given by the definition 5.4 
of [33] and verify on the initial data po £ L^{^) n L^'^°'{^). By considering pon = 
max(po + ^,n) which is strictly positive with boundary condition such that Pn{t, ') = ^ 
on the boundary of O, by the theory of the quasilinear equations (see Ladyzhenskaya 
et al [26] or Friedman in [13J) it exists global classical solutions which verify the energy 
estimates: 

r f \vp^\'dxdt +^ [ p^+'dxdt <^ [ p?+^ 

Jo Jn en -\- i Jq a -\- i Jq 

Furthermore by applying the maximum principle we know that (pn)neN is a decreas- 
ing sequel so to converges everywhere to a limit p. Since pn is uniform bounded in 
L^(L^+"(ri)) for any T > 0, it implies that up to a subsequence pn converge weakly to 
p in L^ {L^~^°' (Q)) . Furthermore Vp" is uniformly bounded in L^(L^(i7)) for any T > 
it implies that up to a subsequence Vp" converges to a limit ip. 

By applying the lemma[T]to p" and the fact that p" is uniformly bounded in L^ {L^^"" (Q,)) 
with e > and converges almost everywhere to p" it shows that p" converges strongly 
to p" in Lj^^{{0,T) X Q,). In particular we deduce that ^p = Vp° 



oCf 



L^ global weak solution on a bounded domain Q, 

By the fundamental L^ contraction principle which ensures that for any solution pi, p2 
of the porous media in L^ we have: 

\\pi{t) - p2{t)\\L^n) < IIpoi - Po2\\LHn)- (2-30) 

The L^ limit solution p of the equation (J1.12p consists in considering the L^ limit of an 
approximate energy solution p„ of ()1.12p with pon converging to po in L^{Q). By the 
estimate ()2.30p we check easily that the limit p does not depend on the choice of the 
regularizing sequel pon- In addition we can verify that p is in C([0, +oo),L^(ri)) (see 
pl29 in [33]). However an important question remains, is the limit solution p is a weak 
solution according the definition 5.2 of [33]. The answer is positive. Indeed for any 
approximative solution we can check that: 

Pn- Pm\{t,x)S.{x)dx+ / / \\p"{s,x) - p'^{s,x)\dxds < / \pon- POm\{x)^{x)dx, 

Q Jo Jq Jn 

(2.31) 
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where ^ is the unique solution of the problem: 

A^ = -1 in [7, ^ = on dQ. 

It implies that p" is a Cauchy sequence in L^((0,T) x $7) for any T > and then 
converge strongly to tp but pn converge strongly in L^{L^{Q) for any T > then up 
to a subsequence r/io„ converges almost everywhere to p and p^ up to a subsequence 
converges to p" and ip then ifj = p°^ which implies that p is a very weak solution of p.l2p 
in the sense of the definition 5.2 of 1331. 



L^ solution on 

We start by approximating the initial data by setting: 

PQn{x) = m.&x.{-n,TQ.\u{pQ{x),n)x{nx), 

whit X £ Co°(M^) such that the support of x is embedded in the ball -6(0, 2) and x = 1 
on -6(0, 1). We have seen that it exists global L} solutions for the homogeneous Cauchy- 
Dirichlet problem in bounded domain ri„ = -6(0, 2n). It suffices to pass to the limit when 
n goes to infinity by using the same type of compactness argument. 
Let us mention to finish that the uniqueness of the L^ solution is a direct consequence 
of the L^ contraction principle. D 

Remark 22 Let us recall that there exists global weak solution which are not classical 
it means not C°° even if the initial data is C°° (see a example due to Aronson in the 
problem 5.7 of JM/)- 

We are now to recall the so called L^ — L°° smoothing effect (as for the dispersive 
equations), we refer to [53] p 202. 

Theorem 2.5 For every t > we have: 

p{t,x) < C||/)o||ii(iRiV)t"'^, 

with a = ivTTTTTo; /^ = j\i( -11+2 ^^^ C > depends only on a and N. The exponents 
are sharp. 

Let us finnish by giving a more general theorem of existence of global strong solution for 
(|2.28p with some properties on the solutions (see [33] p 204-205). 

Theorem 2.6 For every po ^ L"'^(M ) there exists a unique global strong solution of 
(MW such that p G C([0, +oo), L^{R'^)) H L°°{{t, +oo) x M^) with r > 0. Furthermore 
we have the following L°° estimate: 

|p(t,x)|<C||po|IIi(MiV)t-^ 

with a = jv(a-i)+2 ' f^ ~ N(a-i)+2 "^^ C > depends only on a and N. Moreover we 
have the following properties: 

1. The solutions are continuous functions of {t,x) in Q with a uniform modulus of 
continuity for t > r > 0. 
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2. The maximum principles holds. 

3. if Pq is strictly positive and continuous, then p E C°°(Q) Pi C{Q) and is a classical 
solution of li2.28\) . 

4. For every p G (l,+oo) we have the following regularizing effect, p{t,-) belongs in 
LP(M^) and: 

||p(t)||LP(RiV)<C^i~"1|P0|l2?(KiV), 

with ap = (^N^a-l)+2)p ^^d "P = (Ar(„_i)+2)p- 

Let us conclude this section with two important theorem showing the finite speed of 
propagation for the porous media equation (see [MJ p 210) and the time asymptotic 
behavior of the solution which converges to Barrenblatt solutions (see [35] P 69). 

Proposition 2.1 Let p he the global strong solution of \2. 28\) with initial data po G 
L'^{R^) n L°^(M^) and assume that po has a compact support then for every t > the 
support of p{t,-) is a bounded set. 

Theorem 2.7 Let p{t,x) be the unique global strong solution of h2. 28\} with initial data 
po G L^{^ ), pq > 0. Let Um be the Barrenblatt with the same mass as pQ. Then we 
have: 

hm ||p(t)-F„(t)||ii(K^) = 0. (2.32) 

Convergence holds also in L°° norm: 

hm tP\\p{t) - F^(t)koo(KiV) = 0, (2.33) 

with 13= jv(»-l)+2 - 

Remark 23 For more results in this direction we refer also to U^ . 

Let us conclude this section by giving general results (essentially extracted form the 
chapter 9 from [33j ) on porous media equation of the form: 



(2.34) 



f dtp - 2Apip) = 0, 
[p{0,-)=Po- 

2.1.1 General viscosity coefficients 

We assume here that p verifies the following assumptions: 

• /i is a continuous and increasing function: M — t- M with ^(0) = 

• p has at least linear growth at infinity in the sense that it exists c > such that for 
large |s| we have: 

\p{s)\ >c\s\ >0. 

Definition 2.5 A locally integrable function p defined in Qt is said to be a weak solution 
of the problem if: 

18 



1. fi{p)eL'iO,T,H\R^)) 

2. p satisfies the identity: 

/ / (Vfiip) -Vf — pdt<p)dxdt = / po{x)ip{0,x)dx, (2.35) 

J Jqt Jr^ 

for any function tp £ C^{Qt) which vanishes for t = T and has uniformly hounded support 
in the space variable. 

We define L^(M^) by the set of measurable function p^ sucfi tliat p{po) G L^(IR^). We 
sliall consider ^ the primitive of //: 

V'(s) = / p{T)dT. 

Jo 

Let Xt = Lf,{R^) n Li(IR^) and Y = L'^iQr) n L^{Qt) for < T < +oo. 

Theorem 2.8 Let po £ ^ ■ Then it exists a unique global weak solution defined in 
(0,+oo) and Vp{p) G L?{Qt)- Moreover we also have p € L°°((0, r),X). 

2.2 Fast diffusion equations, < a < 1 

The situation is different in the case of fast diffusion equation < a < 1, indeed in 
this case infinite propagation occurs and solution may even vanish in finite time. Let 
us mention that when a is in the interval (0, ?7ic) with ttt-c = max(0, -^—) then it can 
appears a phenomena of extinction of the solution in finite time. In particular it implies 
a lost of the initial mass when po is in L^ (it implies also a lost of the regularity of 
the solution). We refer to |34| theorem 5.7 for a necessary condition of extinction, in 
particular the solution belongs in an appropriate Marcinkewitz space Mp*{M. ). 
Let us mention that in the case a S {rric, 1) the situation is quite similar to the case of the 
porous media equation (except the infinite propagation speed) as the mass is preserved 
which implies no extinction in finite time. Moreover we have self similar solutions also 
discovered by Barrenblatt that we can write under the following form: 

-1 

U^{t,x)=t-^^F{xt-^) with F{x) = {C + Ki\x\'^)l~\ 

with Ki = 2Na ■ ^^ particular the proof of the most of the previous result in the 
last section are based on the existence of Barrenblatt solutions and on the maximum 
principle or in other word the L^ contraction principle. This two fundamental point arise 
also in the case a E {rric, 1) which implies that we have the most of the result of the case 
a > 1 exists also in this case (using essentially the same proof). In particular similarly 
to the case a > 1 in the situation Uc < oi < 1 the global strong solution converges 
asymptotically to a Barrenblatt solution and we have also regularizing effect L^ — L°° . 
For more details in this situation we refer the reader to the excellent books of Vazquez 
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3 Entropy inequality for the quasi-solutions and basic tools 

3.1 Entropy for the quasi-solution of the system ( 11.101) 

We now want to establish new entropy inequalities for system (jl.lOp by applying the 
entropy inequalities discovered in P, ^29] . More precisely if we assume that {p, u) are 
classical solutions of system (jl.lOp . we obtain the following proposition. 

Proposition 3.2 Asswne that {p,u) are classical solutions of system m.lO\) then for all 
t > we have the two following entropy: 



/ -p\u\^{t,x)dx+ / / 2p{p)\Du\^dxdt+ / / X{p)\diYu\^ dxdt 

2 



(3.36) 



Pol^ol {x)dx. 



-p\u + V(p{p)\'^{t,x)dx + - / p{p)\Vu-^Vu\^dxdt 

2 2 Jq JjrJV 



-po\uo + \/ip{po)\^{x)dx. 



(3.37) 



Proof: In order to obtain ()3.36p it suffices to multiplying the momentum equation by u 
and integrating over (0,t) x M^. 

Let us briefly recall the proof of the second entropy (|3.37p introduced by Bresch and 
Desjardins in [6]. To this purpose, we have to study: 

7 /" -I -1 

^ J [^pW\^ + pu ■ Vvj(p) + -p\V^{p)\^]dx. 

step 1: 

First by the mass equation, we have: 

pVuVip{p) ^V(p{p)dx + / p if {p)/S.ip{p)d\Yudx + I p\V^p{p)\ dwudx. 



Step 2 

It remains to estimate the derivative of the cross product: 

— / /)u • Vip{p)dx = / Vip{p) ■ -^{pu) + pu- —Vip{p)dx 

Vv3(p) • -^(.pu) - / div(/)u)v3'(p)— pdx 
V'^(p) • -^(.pu) + / dN{pufip'{p)dx. 
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Multiplying the momentum equation by 'V<p{p), we get: 

Vifip) ■ -f^P'^^ = ~ / (^mC/') + A(p))Av?(p)div'udx + 2 / Vn : Vip{p) ® Vfi{p)dx 



ix, 



— 2 I V(p{p) ■ V ij,{p)divudx — j 'Vip{p)div{pu (g) u)dx, 

where we use the fact that: 

/ V(A(p)divn) • Vip{p)dx = — j X{p)Aip{p)di'vudx, 
and: 

/ div {2 p{p)D{u)) ■Vip{p)dx = / dj{p{p)djUi)diip{p)dx + / dj{p{p)diUj)dnp{p)da 

di{p{p)djUi)djip{p)dx + / dj{p{p)diUj)dnp{p)dx, 
dip{p)djUi)djip{p)dx — / dj p{p)diUidjip{p)dx 
p{p)diUidjjip{p)dx + / djp{p)diUjdnp{p)dx - j dip{p)djUjdnp{p)dx 

p{p)djUjdiiip{p)dx 
2 / Vu : y p{p) ® 'Vip{p)dx — 2 / V/u(p) • 'Vip{p)divudx — 2 / p{p)Aip{p)divudx. 



Step 4 

Sincec/?, p and A satisfies (|1.5|) and (|1.6|) . then we obtain: 

T /" I V7 / \ 1 2 f r 

— ( / pu -Vip^p) + p dx) = — \/ip{p)d\Y{pu^ pu)dx -\- / if {p){d\v{pu)) dx. 

Finally we have: 

Vip{p)di-v{pu0 pu)dx + / If {p){di-v{pu)) dx 

ip {p)u ■ Vpdi-v{pu) — If {p)Vp{pu ■ Vu) + (p {p){div{pu)) dx 
pip (p)divn div (yOu) — pp {p)Vp{u ■ Vu)dx 
p p (p)(divu) + pp {p)u ■ Vpdivti — pp {p)Vp{u ■ Vu)dx, 
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then by ()1.5p and ()1.6p . we get: 

— / V(/j(p)div(yon /9u)dx + / (p {p){di-v{pu)) dx, 

= 2 ph'{p){dwuf + Vh{p) ■ udivu - \7{h{p)){u ■ \7u)dx, 

= 2 pp, {p){divu) — p{p){diYu) — p{p)u ■ Vdwudx + 2 p{p)diUjdjUi 

+ p{p)u ■ VdYvudx, 
{2pp (p) — 2p{p)){divu)'^ + 2p{p)diUjdjUidx 

A(/9)(divn) + I 2p{p)diUjdjUidx, 

whichgives: 

= / X{p){dwu) dx + 2p{p)diUjdjUidx, 
Adding this equahty and (|3.36|) . and using the fact that: 

2p{p)\iu)\^- [2p{p)d,Ujd,u,dx= f2p(p)(^i^^i^^f, 



we easily get (j3.37p . D As in [29] we are also interested in getting a gain of integrability 
on the velocity. We have then the following proposition. 

Proposition 3.3 Asswne that: 

2p{p) + NX{p) > u\{p) 

for some v G (0, 1) (which is a part of lil.6\) ). Then it exists C > such that smooth 
solutions of /11.10\) satisfy the following inequality: 

7 /* 1 I I I ^ /* /* 

J- / P J^ln(l + \u\'^)dx + ^ / p{p)[l + ln(l + \u\'^)\Du\'^dx <C p{p)\Vu\'^dx. 

(3.38) 
for any 5 G (0,2), and with \Vu\'^ = J2iJ2j {diUjl"^ . 

Proof: Multiplying the momentum equation by (l + ln(l + |up))u, we get: 

|p^[^^^ln(l + |n|2)]dx + y"pn.V(i^^ln(l + |n|2))(ix 

+ [ 2pip){l + ln{l + \u\'^)\D{u)\'^dx + f 2p{p)^^^^djUkDij{u)dx 
J J 1 + \u\'^ 

/f 2uuk 

A(/3)(l + ln(l + |n|^)divu|^(ix + / A(p)- — ^j— ^^jUfcdivtidx = 0. 
J i + \u\ 
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Since: 

Idivup < N\Vu\^ and uii{p) < 2/x(p) + NX{p), 

we obtain: 

/'p^[i^^ln(l + |np)](ix+ /'pn.V(i^^ln(l + |np))(ix 



+ z/ /"^(/9)(l + ln(l + |n|^)|L»(n)pdx<C7 f fi{p)\Vu\'^dx. 



(3.39) 



Moreover multiplying the mass equation by — ^^ln(l + |u|^) and integrating by parts, 
we have: 



f i±^ln(l + \u\^)-^pdx - ! pu- V(i^^ln(l + \u\^))dx = 



We deduce that: 



d f l + lnl' 
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-ln(l + |n|2)dx + ^ / /i(p)[l + ln(l + \u\'^)]\D{u)\'^dx 



dtj "^ 2 

<C I ii{p)\Vu\^dx. 

It concludes the proof. D 

3.2 Entropy for the compressible Navier-Stokes equations 

We are going now to consider the following system with p{p) and \{p) verifying (|1.5p : 

' dtp + div(pu) = 0, 

< (9t(pn)+div(pnOu)-div(2^(p)L>u) - V(A(/9)Vn) + eVP(p) = 0, (3.40) 

, (/5,«)(0,-) = (Pc^o)- 

Here P corresponds to the pressure and we shall consider a 7 law P{p) = ap^ with 7 > 1 
and a > 0. 

Proposition 3.4 Assume that {p, u) are classical solutions of system \3.4C^ then for all 
t > we have the two following entropy: 



I \-p\u\'^{t,x) + ^%p'^it,x)]dx+ I I 2p{p)\Du\^dxdt 

JRN 2 7 - i Jq Jf,N 

+ X{p)\divu\'^dxdt = / rpo|^io|^(2;) H -p2{x)]dx. 

Jo Jrn Jf,N 7 - i 



(3.41) 



\-p\u + Vip{p)\'^{t,x) + -^p''{t,x)]dx + - / p{p)\Vu-^Vu\'^dxdt 
2 7-1 2 ./n ./njiv 



+ ea / Vp'^ ■V(f{p){s,x)dsdx = [-po\uo + V(f{po)\'^{x) 

Jo Jrn Jirjv z 

jPoit,x)] dx. 



10 Jm.^ 

ea 



(3.42) 



23 



Proof: We refer to pj for the proof. D 

Proposition 3.5 Assume that: 

2fi{p) + NX{p) > ui\{p) 

for some u £ (0, 1) (which is a part of /ll.6\} ). Then it exists C > such that smooth 
solutions of \1.10\) satisfy the following inequality: 



f P ^^^ ln(l + \u\'^)dx + ^ f pip)[l + ln(l + \u\'^)\Du\'^dx < 



d_ f 1 + W^ 
lit 



2, ^ x^ . (3.43) 



C I pip)\Vu\^dx + C€^ I (^^— — )2-Mx) ( I ip\u\^ + p)dx) 



2 



for any 5 G (0, 2). 

Proof: The proof follows the same lines than the lemma 3.2 in [29], for the sake of com- 
pleteness let us adapt this proof to our situation. Multiplying the momentum equation 
by (1 + ln(l + |np))n, we get as in (fOO]) : 

/ '^di^ ^ln(l + |n|2)]dx+ / pu-V{ j-^ln{l + \u\^))dx 

+ vi I p{p){l + \n{l + \u\^)\D{u)\'^dx< -ae f [I + ln{l + \u\'^)]u ■ V p^ dx (3.44) 

+ C p{p)\\7u\'^dx. 
It remains to bound the right hand side. We have: 
|e f[l + ln{l + \u\'^)]u-Vp'^dx\ 
^ ^1 / -r-prdiUkP^ dx\ + el / [1 + ln(l + \u?')](i\Nup"' dx\, 

J 1 + pr J 

< 2e{ p{p)\Vu\^dx)^ -^dx)^ + e| [I + ln{l + \u\^)]dwup'''dx\. 
Let us deal with the last term on the right hand side: 
e| / [1 + ln(l + \u\'^)]dwup'^dx\ < 

< e( / [1 + ln(l + \u\'^)]p{p){diYufdx)^ / [1 + ln(l + \u\'^)]-^dx)^ , 

We deduce that it exists C > such that: 
|e /"[I + ln(l + \u\'^)]u ■ Vp^dxl < f p{p)\Vu\^dx + y( f [1 + ln(l + \u\'^)]p{p){divufdx) 

,27 



+ |_(|p + Mi + H^)]i^^) 
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where the last term satisfies {if 5 G (0, 2)) for a C > 0: 



927 , ./ r ,n27-i^ 2 \ 2^ 



e^ / ^l^l'^^^-'M/ (^Tj)^^^) " ( / /5[2 + ln(l + |n|2)]ldx)t, 



,27- ' ^ — 



^^''{lij(^)'^dxy-\l{p\u\' + p)dxf^ 



and the proposition follows. D 

3.3 Basic results of compactness 

We would like to finish this section by giving very useful propositions of compactness 
that we shall often apply. We are going to recall the so-called Aubin-Lions theorem. 

Proposition 3.6 Let X ^-^"^ B ^^ Y be Banach spaces (with X which is compactly 
imbedded in B) and (/n)neN 0, sequence bounded in L'^{{0,T), B) n L^{{0,T),X) (with 
I < q < +00J and (^/n)neN bounded in L^{{0,T),Y).Then (/n)neN is relatively compact 
in LP{{0, T), B) for any I < p < q. 

Let us recall now the theorem of Arzela-Ascoli. 

Proposition 3.7 Let B and X Banach spaces such that B M-M- X is compact. Let 
/n be a sequence of functions I ^ B (with I an interval) uniformly bounded in B and 
uniformly continuous in X. Then there exists f £ C^{I, B) such that fn^'f strongly in 
f G C^{I,X) up to a subsequence. 

Lemma 1 Let K a compact subset of M^ (with N > 1) and v"" a sequel such that: 

• v"" is uniformly bounded in L^^°^{K) with a > 0, 

• v'' converge almost everywhere to v, 

then v'' converges strongly to v in L^{K) with v G L^^°^{K). 

Proof: First by the Fatou lemma v is in L^^°^{K). Next we have for any M > 0: 

v"^ — v\dx < I \v'' — v\dx + I \v'^ — v\dx. (3.45) 

K JKr\{\v^-v\<M} J Kr\{\v^-v\>M} 

We are dealing with the second member of the right hand side, by Holder inequality and 
Tchebychev lemma we have for a C > 0: 

/ \v^-v\dx<{i \v'' -v\^'^'^dx)^{{\v^ -v\>M})^, 

JKn{\v^-v\>M} Jk .g^^gX 

< — ^. 



Mi+° 



In particular we have shown the strong convergence of v*^ to f , indeed from the inequality 
(j3.45p it suffices to use the Lebesgue theorem for the first term on the right hand side 
and the estimate p.46p with M going to +00. D 
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Lemma 2 Let f € H" with s > and f e L^ + L"^ with 1 < p < 2. Then f e L'^ . 
Proof: Indeed we have as / G H^: 






so /l{|j^>i} G L2(R^). And as / = /i + /2 with /i e LP(M^) and /a G L^. By using 
the Riesz-Thorin theorem, we know that /i E L'^(M^) with ^ + ^ = 1. As(7>2we then 
have /ln?j<ii £ L^(M ). This concludes the proof. D 

4 Proof of theorem 11.11 

Let us assume in a first time that the solution (p, u) of system (jl.lOp are classical, we are 
going to search solution under the form: (p, — V(/?(p)). The mass equation give us: 

dtp-d\w{p\I^{p)) = Q (4.47) 

2m (P) 

9ip - 2A^(p) = (4.48) 



Since if (p) = ^ ^^' we get: 



Let us check that the second equation is compatible with the first and keep an irrotational 
structure. First we have: 



(4.49) 



dt{pu) = -dt{pV^{p)) = -2Vdt{p{p)). 

-2diY{p{p)Du) = 2div(/i(p)VVv3(p)), 

= 2p{p)VA^{p) + 2Vp{p) ■ VV^{p). 

-V(A(p)divn) = -ViXip)Aipip)). 

Next we have: 

div{pu (g) u)j = ^ di{puiUj) = ^ di{pdnp{p)djip{p)), 

i i 

= ^di{pip' {p)djpd^ip{p)) = 2^di{djp{p)diip{p)), 

i i 

= 2Aip{p)djp{p) + 2{Vp{p) ■ VVip{p))j. 

We have then: 

div(pu 0u) = 2Aip{p)Vp{p) + 2Vp(p) • VVip{p). (4.50) 

Combining (j4.49p and (|4.50p we obtain: 

div{pu (g)u) - 2div{p{p)Du) = 2Aip{p)Vp{p) + 2Vp{p) ■ Wip{p) 

+ 2p{p)VA^{p) + 2Vp{p) ■ VV^ip), (4.51) 
= 2V{p{p)Aip{p)) + 2V(Vp(p) • Vifip)). 



26 



Finally using ()4.49p . (|4.5ip and the fact that X{p) + 2fi{p) = 2/u'(p), we obtain: 
dt{pu) + d\v{pu n) — 2div{p{p)Du) — V(A(p)divn) = 

- V{2dtpip) - 2pip)Aipip) - 2Vp{p) ■ V^ip) - A(p)Avp(p)), (4.52) 

- V{2dtpip) - 2pp{p)Aip{p) - 2Vp{p) ■ Vifip)). 

Next we have since p (p) = \pip (p): 

i i 

= i^Ydi{pdi<f{p)), 

i 

= -pAip{p) + -Vp- Vifip). 

and: 

p'p)Ap{p) = -pAifip) + -p{p)Vp ■ V(p{p), 

= :^P^v{p) + :^^p{p)-^v{p)- 



(4.53) 



In particular from ()4.53p we have: 

Ap (p) = 2pA^{p) + 2Vp{p) ■ Vifip) (4.54) 

Combining ()4.55p and ()4.54p we have: 



dt{pu) + div{pu (g) u) — 2div{p{p)Du) — V{X{p)divu) 

= -V{2p{p){dtp)-2Ap{p))). 



(4.55) 



This concludes the proof inasmuch as via the above equation the momentum equation is 
compatible to the mass equation and must verify the equation (|4.48p . 
But when we take initial density in L} non negative and continuous, we know via the 
theorem 12.61 that the unique global solution of (|4.48p is classical and non negative. It 
justify in particular all the previous formal calculus and prove that {p,u = —Vif{p)) is a 
classical solution of (jl.lOp with p verifying ()4.48p . It concludes the proof. 
Furthermore the different properties on p are a direct consequence of theorem 12.41 12.61 
and [2^ 



5 Proof of the theorems 11.21 

We now present the proof of theorem 11.11 extending to general viscosity coefficient the 
results of p9] in the case of the quasi solutions. Let us begin with recalling the assump- 
tions on the initial data. Indeed we assume that it exits a sequence {pn,Un) of regular 
global weak solution verifying the system (jl.24p (or at least an approximated system, 
typically by using Friedrich approximations). 
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Initial data: 

In particular the initial data Pq,Uq) must uniformly in n satisfy ()1.25p . and ()1.26|) in 
order to verify the entropy inequalities from section [3l more precisely we shall have: 

• p1^ is bounded in L^(R^), p[J > a.e in M^, 

• pJJKI^ is bounded in ^^(R^), 



• 



^Vifip^) = Vfip'') is bounded in L 



2tTn>N\ 



• /)^|'u^pln(l + |ti[J|2) is bounded in L^{R'^). 

With those assumptions, and using the entropy inequalities (j3.36p . (j3.37p and the mass 
equation, we have the following bounds: 



(5.56) 



ll\//'n|lL°°((0,T),L2(]RiV)) < C, 
\\VPn'^n\\L°°({0,T),L^{R^)) < C, 
l|V/(Pn)||L°°((0,T),L2(RAr)) < C, 
\\VPn^Un\\L2({0,T)xR^) < C, 

and 

llPnl^nl ln(l + |Un| ) ||loo((o,T),L1(K^)) < C*- (5.57) 



Remark 24 Let us point out that the gain of integrahility on Un in (5.51) will he a direct 
consequence of a gain of integrahility on the pressure with some restriction on 7, ui and 

V2- 

The proof of theorem 11.11 will be derived in three steps and follows some arguments 
developed in [29j. In the first step, we deal with the strong convergence of the density 
(/9n)neN which enables us to show the convergence almost everywhere of (/3n)neN us to a 
subsequence. We shall also prove the strong convergence of a momentum sequel of the 
form ^ffh\h{^pn)Un with a function h to precise to a function ^h{p)u. In the second step 
we derive the strong convergence of y/]hiUn to y^n in L^g^((0, T) x M) (it allows us to 
give sense to the momentum product PnUn®'^n) by taking advantage of the uniform gain 
of integrahility on un via the entropy inequality ()3.38p . Indeed it will suffice to use the 
lemma[T]after proving almost everywhere convergence via Sobolev injection. In this part, 
we also shall deal with the strong convergence in the distribution sense of the product 
\fPn\fPnUn- lu the last step we will treat the diffusion term which will achieve the proof 
of the theorem 11.21 



Step 1: Convergence almost everywhere on p„ and pnUn 

We are going to begin with proving a technical lemma giving uniforms estimates on p„ 
via the entropy ()3.36p and ()3.37p . 
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Lemma 3 When N = 2,3 V(^^) is uniformly bounded in L^{L'^{R^)) for anyT >0. 



• When N = 3 it implies that ^^, ^^ are uniformly bounded in L°^{{0,T),L^{R^)) 
for any T > which gives: 

pl'^^ is uniformly bounded in L°°(0,r;L^(M^)). (5.58) 

• When N = 2 we distinguish two cases: 

• 1^2 ^ 2, pn is uniformly bounded in L^{L'^{M. )) for any q G [l,+oo[ and any 
T > 0. It implies that ^^ and ^^ are uniformly bounded in L'^{{Q,T),L'i{K)) 



for any compact K . 

• < V2 < 2, Pn is bounded in L^(L'^(M^)) for any T > and any q G [l,+oo[. 
It implies that ^^plJ and -^£sl are uniformly bounded in L°°{{0,T), L'^{K)) for any 
compact K . 

Proof: When A^ = 3, we observe that: 

V(^) = 2M(P)VVP-^VP, 

Vp 2p-2 

= 5V/(rt - ^v,. 

and from conditions ()1.5p . ()1.6p and the fact that /x(p) > we obtain: 

2p{p)p = A(p) + 2p{p) = '^^^^ + '^^^^ > ^Mp). (5.59) 



We deduce that: 



It provides that 



I^Vpl < 2|fiM||Vp|, 
2^02 ^ Vp 

<|;|V/(p)|. 



\v{!^)\<c\vf{p)\. 



It implies by energy estimates that: 



l^( , )llL°°(0,r;L2(RiV)) < C. 

y Pn 



Sobolev's embedding ensured that ^^^ is bounded in L°^(0,T; L®(E^)). Next by using 



!|1.9p we have: 

VPn 
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1 1 I "1 
It implies in particular pn ^ '^^ is uniformly bounded in L^{0,T; L^(]R )) (this is due 

to the fact that pnl{p,,<i} is uniformly bounded in L°°(0,T; L^(R'^))nL°°{0,T; L°°(M^)) 
and that 3 - f + ^ > 1 when N = 3). 



When A^ = 2 by following the same proof than in the case iV = 3 we obtained that 
y^/fi^-j ^g uniformly bounded in L°°(0,T;I 

way (|5.59p when N = 2 and ()1.9p we have: 



y^/fi^-j ^g uniformly bounded in L°°(0,T;L^(M^)) for any T > 0. By using in a similar 



Cpn^^'^'' \Vpn\ < 2|^^^||Vp„| = \Vf{pn)\ when Pn > 1, 

\/Pn 

Cpn'^'^lVpnl < 2|^^^||Vp„| = \Vf{pn)\ when Pn < 1. 



(5.60) 



When ui > 2, choosing ip G Cq°(M^) with ^ = 1 on -6(0, 1) and supp^ included in B{0, 2) 
we have: (1 — il^ipn))"^ y/Pn is uniformly bounded in LI^(L^(M^)) for any T > 0. Since 
y^ is uniformly bounded in L!^(L^(M^)) for any T > we deduce that (1 — 1p{pn))^/Pn 
is uniformly bounded in L^(ii'^(R^)) for any T > 0. It implies that (1 — ilj{pn))pn is 
bounded in L^(L''(R )) for any q G [l,+oo[ by Sobolev embedding. Let us deal now 
with the term ipipnjPn which is bounded in L^{L^{R^)) D L^(L°°(IR^)), it proves that 
Pn is bounded in Lf{Li{R^)) for any q G [l,+oo[ . Via ([HD It implies that ^^ and 

^^' are uniformly bounded in L°°((0,T),L''(M^)) for any compact K. 



Let us deal with the case < z^2 < 2. By using ()5.60p we show that V((l — il^{pn))pn ) 
is bounded in L^{L'^{R^)) and (1 - i;{pn))pt is bounded in L~(L^(M^)). Since by 

ill 

Tchebytchev lemma (l—ip{pn))pn is strictly positive only on a set of finite measure it im- 

plies that {l-7p{pn))pn is also bounded in L^{L'^{R^)). We deduce that (l-V'(pn))/On 

ill 
is bounded in L^(if^(]R^)). By Sobolev embedding it yields that (1 — 'ip{pn))Pn is 

bounded in L^(L''(R^)) for any T > and any q G [1,+cxd[. Since ip{pn)Pn is bounded 

in L^(L~(M^)nLi(M^)) for any T > we conclude that pn is bounded in L^(L«(M^)) 

for any T > and any q G [1, +cxo[. 

The proof is similar for ^£— ^ by using the remarks [71 D 

Lemma 1 If p{p), X{p) satisfies lll.5\) . Ill.6\) and in addition we assume that g(x) = ^-^ 
is a bijective function on [0, +oo) and that g~^ is continuous, then when we distinguish 
the two following cases, we have: 

• i/i > 2 

1. !^ is uniformly bounded in L°°{0,T;HI^{M.^)) 

2. dt>^ is bounded in L'^{Q,T-W[^l''^{W^)). 



Pn 

As a consequence up to a subsequence (via the Cantor's diagonal process) i^-^^J- 
converges almost everywhere and strongly in C([0,T], Lf^^{R ))) to v. In particular 
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it implies that ^^^' converges up to a subsequence almost everywhere to v and we 
define p as follows: 

It implies that pn converge up to a subsequence almost everywhere to p. 

• 0<ui <2 

Let us consider: 

with < «! < mm(i, j^)- 

1. I3{pn) is bounded in L'^{0,T; H^{R^)) 

2. dtP{pn) is bounded in L^{0,T;W-^'^{R^)). 

As a consequence up to a subsequence, /3{pn) converge almost everywhere and 
strongly in C{[0,T],Lf^^{R^))) to v. We define p by: 

p = r\p). 

In particular we have that pn converge up to a subsequence almost everywhere to p. 

Furthermore pn converge strongly to p in C([0,T], Lj^"(M^)) if N = 3 with a > small 
enough and in C{[0,T],L^^^(W )) for any p > 1 if N = 2. This last result is under the 
following hypothesis: 

• When 2 + N < i^i, we assume that g and g are increasing on (0, +oo). 

Remark 25 Let us point out that we could weaken the last assumption on g when z^i > 2 
by assuming that g and g are only increasing in a neighbor ofO and 0/+00. As mentioned 
above, let us point out that this last condition is quite natural since this is true when we 
set p{p) = pp'^ with a > |. 

Proof: Let T > 0. We are going to distinguish two case when 1^1 > 2 and when 
< i/i < 2. 

•ui > 2 



The first estimate is a direct consequence of the lemma [3] in the appendix. Indeed 
we know that V ^^ is bounded in L!^(L^(M )). Easily we deduce by lemma [3] that 

is uniformly bounded in L^{L'^{K)) for any compact K. Next we observe that: 



r.tP{Pn). ,. fP{Pn) ■. , f ' f ■. , — 3/i(p„)s 

y/Pn y/Pn 2 y/pn 

I (5.61) 

= -div( y/p^Un) + ( , ,\ o y ) \//"(Pn)dlVn„. 

Pn VPiPn) ^^ Pn ' 
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Let us start with estimating the first term on the right hand side of ()5.6ip . According to 
we know that: 



"^^^"^^^^"^ < C^{pn)pr~^ when pn < 1. 

Pn 

Since ^ - ^ > it imphes in particular that ^(P")^(P") is bounded in L^(L°°(M^)). In 
order to have local estimates, it suffices to deal with the region when />„ > 1, it means 

(l-V>(pn))At(Pn) 
pn 

It suffices to use the lemma[3] which insures that when A^ = 3 '^^pL' belongs in L^(L®(M )) 
and we deduce that (i-V'(Pnj)M(pn) -^ bounded in L^{L^{R'^)). Finally we have obtain 
that ^^1^"' is bounded in L^{Lf^^{M. )) when A^ = 3. It yields the uniform boundness of 
^^^"•^ /]hiUn in LJ^(L^^^^(M^)) when A^ = 3. A similar argument insure the same result 



Pn 

when N = 2. 



By proceeding similarly we also prove that -v / ^^^^^ y7x(/)^divu„ is uniformly bounded in 



pn 



L^{Lj^^{R )). Indeed following the same argument than for the previous term, we know 



that for A^ = 2, 3 via the lemma[3] ^ t!±Ei}d. is bounded in L^(L^'^(K)) for any compact K 



Pn 



and since ^^/ p{pn)dwun is bounded in L|,(L^(]R^)), Holder (s inequality give the desired 
result. 



Let us now deal with the term — / \/ P'{Pn)divUn. A simple calculus using (ll.5p give 

yp{pn) 
us: 

P{Pn)^/p^ _ I, Xjpn) r, IPipn)^ 



VKPnj 2 yjpnp.{pn) V ^" 

We have only to deal with the term . ^^"-^ =, the other one has been estimated. By the 

yPnP{Pn) 

remarks [8] we know that it exists C > such that: 

\\{p)\<Cp{p) Vp>0. 
It implies that: 



Kpn) I ^ ^ lp{Pn) 



\/PnP{Pn) V P"^ 

And as we know that ^^i^ is uniformly bounded in L'^{L^^^{MJ^)) for A^ = 2, 3 via the 
lemma El it achieves the proof of the second estimates. 

By the Ascoli's theorem, the fact that the application u ^ (pu with cj) G C^(M ) is 
compact from H^{M.^) to L^(M^) and the Cantor's diagonal process it entails that ^^p^' 

converges strongly up to a subsequence in C{[0,T],Lf^^{E.^))) to vi = ^^ (and in 
particular in Lf^^{{0,T) x M^) ). We shall define p in the sequel by: 



Vp 



PJPn) 



Furthermore an immediate consequence is that up to a subsequence ^^pL' converge almost 
everywhere to ^-^- And since g~^ is continuous, it implies that up to a subsequence p„ 
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converges almost everywhere to p. 

• < i/i < 2 

In this case we are going to study: 

/3(Pn) = HPn)y/pZ+{l " V'(Pn))p"S 

with Qi to choose suitably. Let us start with dealing with ip{pn)y/Pn, we have then: 

fJ- (Pn) ^JPn 

The first term is easily bounded in L^(L^(IR^)) via the entropy (|3.37|) since the support 
of ■0 is included in the shell C(0, 1, 2). The second term is also bounded in L'^ijJ'iW^y) 

^^^ IpjPn) ;„ Kr^nnrlo^ in TOO l' TOO / TCP ^^ 

Indeed we have 



by observing that 4fH is bounded in Lf{L°°{M.^)) via (US]), ([13]), (US]) and vi < 2. 

M (pn) 



; / J ^ tPiPn Pn 
P [pn) 

The conservation of mass provides that H/Onllii = IIPollLiwhich implies the L°°(0, T;//^) 
bound on ijj{pn)^/fhi- Next we observe that: 

1 / 3 

dt{lp{.Pn)^/K) = -div{'lp{pn)^/p^Un) + (-V'(Pn)V^ " ^ (/On)Pn )divu„. (5.62) 

The first term on the right hand side is obviously bounded in L^(VF~^'^(R )). The two 
last terms in (|5.62|) are bounded in L^(M^~^'^(]R^)). It implies finally that tp{pn)\/Pn 
are uniformly bounded in L^(iJi(M^)) and dt{ilj{pn)^/lhi) in L\{W~^''^{R'^)). 
Let us deal now with the term (1 — ip{pn))Pn^ , we have then: 






V((l - i^{pn))p'^') = -= ^'{pn)p7'^Pn + (1 " ^{pn))\^, ^^Vp„. (5.63) 

P \Pn) JP^ 



By ([rs]), dTj]) and ([TH]) we show that: 



Pn 1^/1 ,;.^. N\ „"i-|+:^-|^ 



1(1 - ?/'(/5n))^7 ^1 < (1 - ^{Pn))Pn 

P [pn 



Since < ai < ^ it implies that (1 - ^(pn))4^ is bounded in LS?(L°^(M^)) and 
(1 - V'(Pn))4^ ^^f^Vpn is bounded in Lf{L'^{R^)). The first term on (f533]> is easy 

/^ \pTi) y Pi^ 



V'(P„) v^ ^^" 
to deal with. Now we have: 



dt{{l - HPn))Pn') = -div((l - ^lJ{pn))pTnn) " ("V' (/>n) 



rn 



+ (ai - 1)(1 - V'(Pn))^F==)\/MPn)divu„. 

\/P[pn) 
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The first term on the right hand side of ()5.64|) is bounded in L^{W ^'^(M. )) since 
(1 — ilj{Pn))Pn^ is bounded in L°°{L'^{M.^)) because ai < ^. We have now according to 

(USD 

1(1 - v^(p„))-£n_| < (1 - ^(p^))^^i-i+2^-^. 

Vf^iPn) 

Since ai < 3 it implies that this term is bounded in LJp(L°^(M )) which shows the 

bound of (1 - tp(pn)) f\ A/;u(pn)divnn in L%{l?{n^)). The second term of EMI is 

easy to treat. 

Finally we have proved that (1 — ilj{pn))Pn^ is uniformly bounded in L^{H^{^ )) and 

dt{{l—ip{pn)) Pn^) is also in L'^{W~^''^ (R'^ )) . In conclusion it shows that /3(p„) is bounded 

in Lf{H'^{R^)) and dt{P{pn)) is bounded L^W'^^^iR^)). 

Thanks to the Ascoli's theorem and the Cantor's diagonal process it gives the strong 

convergence in 

C7([0,T],lL(M^)) 

of /3(/9ra) to V2. It implies in particular up a subsequence of the convergence almost 
everywhere of P{pn) to ^2- We shall define p in this situation by: 

p = r\p), 

indeed /3 is inversible and we verify by continuity of /3~^ that pn converge up to a subse- 
quence almost everywhere to p. 

Strong convergence of pn 

We are now interested in proving the strong convergence of pn to p. Let us deal with the 
first case z^i > 2. 

First case: 1^1 > 2 

For the moment we have only obtained strong convergence on ^^^ to ^-j^, let us trans- 
late this strong convergence on /9„. We are going to distinguish two different cases when 
z^i > 2, let us start with the first one. 

• We assume here that (i^i^) ^ t!-^ are increasing on (0, +00). We have then the following 
lemma. 

Lemma 4 Let gi a regular function with gi(0) = 0. When x,y > and gi, g^ are 

increasing, we have: 

gi{\x-y\)<\gi{x)- gi{y)\. (5.65) 



Proof: It suffices to study the function: 

p{x) = gi{x-y)- (51 (x) - gi{y)), 
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when X > y. We observe that for all x > y > 0: 

p'{x) =g'i{x-y)- g'i{x) < 0. 

It implies that when x > y, p is negative andp is decreasing on [y, +00) with is equivalent 
to say that: 

9i{x -y)- {gi{x) - gi{y)) < p{y) = gi{0) Vx > y. 

It implies that for all x > y we have: 

9i{\x - y\) = gi{x -y)< {gi{x) - gi{y)) = \gi{x) - gi{y)\. 
By proceeding similarly when < x < y we obtain (|5.65p . D 

In particular since we assume that (^^i^) and (^^^) are increasing on (0, +00) we deduce 
from the lemma |U 

K\P- Pn\) ^ Mpn) _ P{p) , 
\/\p- Pn\ ~ \fP^ \fP 

Using p.9p we obtain: 



|2-lv + 2iv < |^1^_^| V|y9-p„| >1 

C ^/p^ ^/p 

2 N + 2N <—\^—J- ^ V/9-p„<l. 

C y/p^ y/p 



(5.66) 



Since y^ converges strongly to ^^ in C([0,T],L^^^) and !^ , ^^ are bounded via 
lemma[3]in LJ^(L^(]R )) when A^ = 3 (it suffices to apply Fatou's lemma), we deduce by 
interpolation that ^-y=- converges strongly to ^^^ in C([0,T],L^J^°) for any small a > 0. 



Similarly when N = 2 we obtain by interpolation and via lemma [3] that ^^^' converges 



/Pi 
strongly to ^ in C{[0,T],Lf^J for any p > 1. 

By using (|5.66|) we deduce that {l—i/j{p—pn))\p—pn\^~^^^ and ^(p— p„)|p— p„| 2^iv + 2iv 
converge strongly to in C([0,T],L;^") for any small a > when A^ = 3 and in 
C([0, r], Lf^J for any p > 1 when N = 2. Since 1^1 > it yields that: 

sup IIIp - />nr^'^^^^-*H{|p-p„|>l}(t, OIIl^W ^n->+oo 0, (5.67) 

te[o,T] 

for any compact K with p(A^) = 6 — a for any small a > and p{N) = p for any 
p £ [1, +oo[ when N = 2. And similarly we have: 

sup IIIp - Pnr(^)(^~^ + ^H||,_^„|>i|(t, ■)\\lHk) ^n->+oo 0. (5.68) 

te[o,T] 



It implies since 1^2 > that pn converges strongly to p in C([0,T],L^^"(M^)) for a > 



small enough when A^ = 3 and pn converges strongly to p in C([0, T], L^^^{M.'^)) for any 



p> 1. 
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• U2<N + 2 

Let us deal now with the case when (^^i^) and (^^i^) are not increasing on (0, +00) 
By calculus and using (jl.Sp we have: 

^ Vp Vp 2 ^t ' 

_ 1 A(p) + /i(/>) 
2 ^ 

Next by (|1.9p we obtain: 



p2 



Cp^-^-^ < |(^)'| < Cp^~^-7r yp< 1, 
C7plw-i-^ < |(^)'| < Cp^-^-TT Vp > 1. 

Vp 



(5.69) 



Assume that ^ — ^ — ;^ ^ which is equivalent to 1^2 < N + 2. In this case we obtain: 

l(^)'l>CVp<l. 

Vp 

Let us recall that the derivative of the inverse function of g{p) = ^-4^ is: 

9 {9 Hp)) 

In particular it means that the inverse function g~^{p) is Lipschitz on the region p < 1 
and more generally on the region p < M, it provides then the following inequality for 
Cm > depending on M (since when p < M we have that ^^ is also bounded via ()1.9p 

and the hypothesis | + ^ — ;^ ^ 0): 

l^=- 7^|l{p<i}u{p„<i}<-Af > \g (— =-) - 9 {^=^)\Hp<M}u{prr<M} 

y/P yPn \JP \/Pn (5.70) 

> |p-Pn|l{p<l}U{p„<l}- 

We deduce the following estimate for any compact K: 



Pnit,x)dx - p{t,x)\dx = / \pn{t,x)dx - p{t,x)\l{pyM}U{p„>M}it,x)dx 

^ ^^ (5.71) 

Pn{t, x) - p{t, x) I l{p<M}U{p„<Af } {t, X)dx 
K 

The second term converges uniformly on (0, T) to when n goes to infinity via ()5.70p 
applied to M > and the strong convergence in C{[^,T],Lf^^{M.^)) of ^^ to ^. Let 
us deal with the first term, since we know via the lemma [3] that p„ is uniformly bounded 
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in L!^(L^^'^(]R )) with e > we have by Holder's inequality and Tchebytchev lemma for 
C>0: 



K 



\pn{t, X)dx - p{t, x) I l{p>M}U{p„>M} (i, X)dx 



< 



( / \pn{t, X)dx - p{t, x)|l+n{p>M}u{p„>Af}(*, X)dx) i^' \{p > 1} U {p„ > 1}| 1+^ 

Jk 

< 2C- 



K 

"poIIli 



M 

This last term goes uniformly on (0, T) to in n when M goes to infinity. It show the 
desired result. 

Second case < z^i < 2 



In this case it suffices to apply exactly the same argument than in the case vi > N + 2 

y/P 



except at the place we consider not ^^^ but 13 (p). It concludes the proof of the lemma. 



D 

Lemma 5 Let ip G C^(M^) with ^ = 1 on B{0, 1) and supptp is included in i?(0, 2). 
We are going distinguish two cases: 

• When vi>2 we set: 

Vn = tp{Pn)pipn)Un + (1 " tp{Pn))PnUn, 

we have: 

• Vn is uniformly bounded in L'^{0,T,W^'^{K) for any compact K. 

• dfVn is uniformly bounded in L'^{0,T,W^'^'^{K) for any compact K. 

Up to a subsequence, the sequel Vn converges strongly in L^(0, T; L^^^(]R )) to some v{t, x) 
for all p G [1, |). In particular: 



(V'(Pn)/w(Pn) + (1 — fp{Pn)))un -^ V almost everywhere (x,t) G (0,T) x 



tN 



Note that we can already define u{t, x) = , . s . \\}i_^( \\ outside the vacuum set 
{p{t,x) = 0}, but we do not know yet whether v{t,x) is zero on the vacuum set (in 
particular if there is no concentration phenomena for v on {p{t, x) = 0}). 

• When < vi < 2 we consider: 

Vn = i^{pn)Pn^'^Un + (1 - ^(/9„))p(^+^U„, 

with (3 < ^ and /3i verifying the following assumptions: 

/3i>l, 

1 ZV2 

^^ + iV - 2iV - °' (5-72) 

1 1 1^2 

we have: 
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• Vn is uniformly bounded in L'^{0,T,W^'^{K) for any compact K. 

• dtVn is uniformly bounded in L'^{0,T,W^'^'^{K) for any compact K. 

Up to a subsequence, the sequel Vn converges strongly in -L^(0, T; L^^^(R^)) to some v{t, x) 
for all p ^ [1,^). In particular: 

{i^{Pn)Pn~^^ + (1 - i){pn))p'^^^)un -^ V almost everywhere {x,t) G (0,r) x M^. 

Note that we can already define u{t, x) = ^ +i^i'i k \\ p+i outside the vacuum set 

{p{t,x) = 0}, but we do not know yet whether v{t,x) is zero on the vacuum set (in 
particular if there is no concentration phenomena for v on {p{t, x) = 0}). 

Proof: Let us start with proving that Vn is uniformly bounded in L^(0, T, VF^'^' 



• We are going to start with the case z/i > 2. 

Let us prove that f^ is bounded in L^{0,T,W^'^{K) for any compact K. We have then: 



i''{Pn) p'jPn) 
1 1 \P^ I 

l -V^(Pn)) P jPn) r, ^- , (1 - V'(/9n ))Pn 

P [Pn) V^ VPiPn 



Ci((l - WKPnjjPnUnj) = 77 TPn -^OiPn^J PnUnj 



yi-ip\Pn)) P {Pn) a ^- . 1,1 - ^P{pn))Pn r^-^r. 
"^ ..//„ N -^l^OiPny/PnUnj H / . . V P{Pn)OiUnj ■ 



We observe that by using (|1.9|) . (|1.5p and i^i > 2 that ^ — JPyP")' jg uniformly bounded in 

M (pn) 

L2?(L°°(M^)). It implies that the second term ^^'/f^^"^^ ^^-^diPn./Kunj is bounded 

-^ M \Pn) yPn ' V ' J 

by Holder's inequality in L|,(L^(R )). The first term on the right hand side of (|5.73p is 
bounded in Lj,(L^(M^)) by using the fact that the support of ip is included in the shell 
C(0, 1,2). The last term is also bounded in L'^{L^ {W^ )) . Indeed by (US]) and z^2 > 2 it 
suffices to observe that: 

(1 - V(Pn))Pn ^ r- r- 



VpiPn) 

which implies that (i-^(P"))P" ig bounded in L^(L^(R^)). 
y/piPn) 

Finally we have seen that Vf^ is bounded in L|,(L^(]R^)). And since we see easily that 

vl is also bounded in Lf{L^{R^)), it implies that v^ is bounded in L'^{W^''^{R^)). 

Let us deal now with estimating dtv'^. A simple calculus gives for any regular function g: 

dt{g{Pn)PnUn) = -g{Pn)PnUn ' "^Un " 9 (pn) pl.UndivUn - UndlV {g (pn) PnUn) 
+ 2div {g{pn)p{Pn)DUn) - 2Vg{pn) ■ p.{pn)DUn + V {g{pn)\{Pn)dlVUn) 

- Vg{pn) X{pn)divun, (5.74) 

= -div {g{pn)PnUn ® Un) - g (pn) plundlVUn + 2div {g{pn)p{Pn)DUn) 

- 2Vg{pn) ■ p{pn)Dun + V(5f(p„)A(/9„)divn„) - V g{pn) A(p„)divn„. 



When we apply the previous formula to g{pn) = (1 — i^{Pn))\ ^ , we have to estimate 



KPn) 



Pn 



all the terms on the right hand side of (|5.74|) . it comes for T > 0: 

g{Pn)PnUn «> ti„ = (1 - ^{pn))\/7hiUn ® ^/p^Un. (5.75) 
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By Holder's inequality we obtain that {'^—•4>{pn))y/PnUn®y/PnUn is bounded in L!^(L^(M^)). 
Next we have: 



3 

2 J- _ I'f \ P^ 



g {pn)plundivun = -ip (Pn) . " - V/^^n\//^(Pn)divn„ (5.76) 

Vf^iPn) 



This term is bounded in L|,(L (M )) by Holder's inequality since i/j is supported in the 



sheh C(0, 1,2) which implies that V {Pn)—r^ — is bounded in LS?(L°^(M^)). Similarly 
we have: 



g{Pn)p{Pn)DUn = (1 - i^iPn))^/ p{Pn)\/ P{Pn)DUn, 

= (1 -ip{Pn))\ PnVp{Pn)DUn 

y \/Pn 

According lemma [3] we deduce easily that (1 — '4'{Pn))\/ p{Pn)\/ p{Pn)Dun is bounded in 
L'^{L^{K)) for any compact K. Next we have: 

'^'^giPn) ■ p{pn)DUn = -2Tp (/9„)Vp„ • p{pn)DUn. (5.78) 

This term is bounded in L|i(L^(M )) by Holder's inequality since ip is supported in the 
sheh C(0, 1,2). The two last term in (fHTD are similar to treat. (fHTTSD . (fHTTUD . (f^TTTD and 
()5.78p implies that 5jf^ is uniformly bounded in L'^{W~^'^{K)) for any compact K. 

• Case < z^i < 2. 

In this case we are going to consider: 

^n = (1 - 'P{Pn))p^PnUn- 

Let us start with proving that Vf„ belongs in L'^{L}{K)) for any compact K. We have 
then: 

«/n If \\ \ '4''{Pn'^^) p'{Pn)r. , 

di((l - 1p{pn))p'nPnUnj) = 77 ^Pn -^^ Oi Pn^/p^Unj 

P- [Pn) \/P^ 

, (1 - V^(/9n))p« p'jpn) r, ^- ^ il-ij{pn))pn'^^ rT^F) 

P (Pn) ,/p^ Vp(Pn) 

By using (jl.6p and ()1.5p . we obtain that: 

, 1 2p{p) + JVA(p) , ,^ 1 Mp) 
^(^^ = 2^ ^ + ('-iv)^' 

>(1-1 + ^)^. 
We deduce that according to (|1.9p that: 



(5.79) 



^ il-lPipn))p^n ^ ^ ^^ (l-^(/9n))p^^ 
p'(Pn) ~ P{Pn) 

(l-V'(Pn))p^\ 



<C|- 



„ N^2N 
Pn 



<C|(l-V(/>n))pn"'^ ^"'^l 
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Since /? < ^ we deduce that (^ ^M)p^ jg uniformly bounded in L^(L°°(M^)). In par- 
ticular it implies that ^^~tl^"!^^" ^^=-diPn./Kunj is uniformly bounded in L^(L^CR^)). 

Since the support of ij: is included in the shell C(0, 1,2) we easily observe that the 
first term on the right hand side of ()5.79p is uniformly bounded in L^{L^{M. )). The 
last term is also bounded in L|,(L^(M^)) since frac{l — ip{pn))Pn \//^(Pn) belongs in 
L^(L2(M^)) by using ([131), /3 < -^ and the lemma [3 

Finally we have seen that Vf^ is bounded in L|,(L-^(M^)). And since we see easily that 
vl is also bounded in Lf{L^(R^)), it implies that vl is bounded in Lf,(P^^'^(R^)). 

Let us deal now with estimating dtv'^. It suffices to deal with the formula ()5.74p and 
replacing g{pn) by (1 — '4>{pn))pn- Let us start with the first term of (|5.74|) : 

g{Pn)PnUn «> n„ = (1 - 'lp{Pn))PnVPnUn '^ \/%Un- (5-80) 

Since (3 < — -^ it implies that {l—ip{pn))pn is L'^ bounded and then by Holder's inequality 
we obtain that (1 — ipipn)) Pn\/PnUn ^ \/7>nUn is bounded in L^{L^{W )). Next we have: 



Prl 



^+13 



^^^^ , (5.81) 

+ /3(1 - ^lJ{pn))—^==^/p^Un^/piPn)dwUn. 

The first term is bounded in L'^{L^(M. )) by Holder's inequality since tp is supported 
in the sheh C(0, 1,2) which implies that ip ipn) ^," is bounded in LS?(L°°(M^)). 



h,8 



The second term is also bounded in L\.{L^(MJ^)) since (1 — tpipn)) ^", is bounded 

yKpn) 

in Lf{L°°{^^)) by using ((L9]) and the fact that /? < ^. 
Similarly we have: 
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{Pn)p{Pn)DUn = (1 - ^{pn))\/ p{Pn)pi\/ p{Pn)DUn, 
= (1 - ^P{Pn))\|^pt^V^4Pn)DUr. ^^'^^^ 

According lemma [3] we deduce easily that (1 — ip{pn))^/ p{Pn)y/ IJ'{Pn)Dun is bounded in 
L'^{L^{K)) for any compact K. The last term gives: 

2Vfi((p„) • p.{pn)DUn = -2lp {pn)p^Vpn ■ p{pn)DUn 

, n,-, II ^^Pn W P(.Pn) P ipn) ^ r-, T^ ^^''^'^' 

+ /3(1 -1p{pn)) rr-. ^Vp„ • y/p{pn)DUn. 

P KPn) y/Pn 

The first term is bounded in Ly(L^(M'^)) by Holder's inequality since ijj is supported 
in the shell C(0, 1,2) and we proceed similarly for the second by observing that (1 — 
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1 



ip{pn) Y" ^,/^^ is L°° bounded via ([LSD, (USD and ([L9]). The two last term in (l57i]l 

M (Pn) 

are similar to treat. ()5.80p . ()5.8ip . ()5.82p and (|5.83p implies that dtv"^ is uniformly 
bounded in L'^{W^^''^(K)) for any compact K. 

Let us deal now with the term v\ . We are going to distinguish two cases. 

• Case vi > 2. 

We have: 

// N / N 1p{Pn)Kpn) ^ 

WKPnjP-KPnj'U.n = ■p= VPnUn, 

VP^ 

where via the lemma [3] fiBjij^Mhil jg uniformly bounded in L^(L°^(M^)) and ^/p^Un is 



uniformly bounded in L°°{0,T;L'^(M.^)) which implies that 'ip(pn)p{Pn)un is bounded in 

L'=°{0,T;L'^{R^)). 

Next we have: 



di{lp{pn)p{pn)Unj) = 1piPn)VpiPn)VpiPn)diUnj + 1p{pn)p Pn)diPnU 



nj, 



^{Pn) VpiPn) VpiPn)diUnj H "^ " diPny/K.Unj. 

V Pn 



i'{Pn)lJ. (Pn) Pt „ \^ „-„\fr^^-^y^rUr^-.-,-^Ar^A\-^ T «= / T 2 I'Trj) TV \ 



By entropy inequality ()3.37p we know that '^^^"'^^"' diPn is uniformly bounded in L^(L (M )) 



which implies that ^^''"'J^^''' ^iPn^/fh^Unj is uniformly bounded in L^(L^(M^)). 



i^nj 



Let us deal now with the term ip{pn) \/ p{Pn) \/ piPn)diUnj , we know that \/p(pn)d-, 
is uniformly bounded in Ly(L^(]R )). Next we know that 'ipipn)^/ p{Pn) is bounded in 
L^(L°^(IR^)) which provides on ijipn)^/ p{pn)^J p{pn)diUnj a L\{L\M.^)) bound. Hence 
for any compact K: 

^{ip{pn)p{Pn)'Un) is bounded in L {0,T;L (K)) 

In particular we have obtained that for all compact K: 

Vj^ = ■ip{pn)p{pn)un is bounded in L {0,T;W ' (K)). 

We are now going to estimate dtv^, it suffices to estimate each term on the right hand 

(Pn 
Pn 



side of iK7^ by replacing g{pn) by V(Pn)^^^^ 



We have then: 

g{Pn)PnUn (^ U^ = V'(Pn) —^/PW^i-n ® ^/pW^n, (5.84) 

Pn 

By Holder's inequality we obtain that ip{pn)p{Pn)un ® Un is bounded in L^{L^{R'^)). 
Indeed we have used the fact that by p.9p and since z^2 > 2 then i'ipn) p" is bounded 

inL^(L°°(R^)). 
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Next we have: 



^^^ (5.85) 

= {tp' {Pn)\/p{Pn)Pn + iJ{pn)^J==)^/P^Uny^p{p^divUn 

The first term is bounded in Lj,{L'^(M.^)) by Holder's inequality since tp is supported 

in the shell C(0, 1,2) which implies that ip (pn) \/ p{Pn) Pn is bounded in L^(L°°(M^)). 
The second term is also bounded in L^(L^(M^)) because using (|1.9|) we observe that 



V'(/Qn)^ / " , is bounded in L^(L°°(M^)). 
The third term gives: 



g{Pn)p{Pn)DUn = 1p{Pn) —^Jp{pn)DUn, 

. ^" (5.86) 

= IpiPn)) Z- \ —y/p{pn)DUn 

xfK. V Pn 



By (dJl) and i/i > 2 we know that v/^^ is bounded in L5?(L°°(M^)), since via lemma[3] 
^^^ is bounded in L^{L^{K)) for any compact K we deduce that ip{pn) p" p{Pn)Dun 
is uniformly bounded in L'^(L^ (K)) for any compact K. 
Next we have: 

'^'^g{Pn) ■ p{pn)DUn = 2tp (/9„)V/9„ • p{pn)DUn 

a(p\ (5.87) 

+ 2lp{pn)y/p{pn)Pn -p^V Pn ' \//U(p„)L>U„. 

\/Pn 

The first term is easily bounded in L2,(L^(M^)) since the support of tp is included in 
C(0,l,2). The second term is also bounded in L|,(L^(R^)) by Holder's inequality be- 
cause via ()1.9p we deduce that ■ip{pn)'\/p{pn)Pn is bounded in L!^(L°^(M^)). The two 
last term in (|5.74p are similar to treat. (|5.84p . (j5.85p . (|5.86p and (|5.87p imphes that dtvj^ 
is uniformly bounded in L'^(W~^'^(K)) for any compact K. 

• Case < i/i <2. 

We are going to work with: 

We have then: 



vi = 1p{Pn)PnPnUn- 
i^{Pn)ptPnUn = 1piPn)Pn ^ y/p^Un- 
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Since /3i > -^, it implies that ip{pn)Pn PnUn is bounded in L!^(L^(M^)). 
Next we have: 



+ MiPn) ," ^ -p^^iPn^/p^Unj + ^j' {pn) Pn^"^^ diPnUnj ■ 

P \Pn) ^/P^ 

Using dm, ([rSD and dTU]) we show that: 

\'^{Pn)^^A<HPn)Pn^^'^ ■ 
P [Pn) 

Since /3i + ^ - ^ > it imphes that V'lpn)-^ is bounded in L^(L°°(M^)). We 

deduce that the second term on the right hand side of (j5.88p is uniformly bounded in 

L^{L^[M. )) . The third term is easy to treat by using that the support of V' is a shell 

C(0,l,2). 

Let us deal with the first term of ()5.88p . By using (jl.Op we get: 

I , — 7- I ^ iP{Pn)Pn 

VPiPn) 



Since /Ji + i + n^ — Tif>Oit provides a L°° bound on ^^P"'P" which insures the 



Lf,(L2(R^) bound of ^(P"W"^' y/^(p„)g^n„^.. Hence for any compact K: 

V(^(Pn)p^'+^^in) is bounded in L^{Q,T;L^{K)) 

In particular we have obtained that for all compact K: 

vl = V'(/>„)p^i+^'u„ is bounded in L^{id,T-W^^^{K)). 

We are now going to estimate dtv\, it suffices to estimate each term on the right hand 
side of (|5.74p by replacing g{pn) by V'(/0n)/0n • 
We have then: 

g{Pn)PnUn ® Un = ij{Pn)pt~^VPn'>J'n <^ ^/P^Un, (5-89) 

By Holder's inequality and the fact that /3i > 1 we obtain that '>p{pn)PnUn Un is 
bounded in L^{L'^{R^)). Next we have: 

3 
„2 

g {pn)plundivUn = {lp {Pn)Pn^ + f3iTp{pn)Pn^~^)—F==\^Uny^p(pn)diyUn, 

^^ (5.90) 

= ii^' iPn)^== + Pli^iPn) !^, - )y/p^Uny^ p{pn)diYUn 

VPypn) vpypn) 
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The first term is bounded in Lj,(L (M )) by Holder's inequality since V is supported 



in the shell C(0,l,2) which show that V {Pn)-^ is bounded in LS?(L°°(]R^)). The 

second term is also bounded in L'^{L^{M. )) because using ()1.9p we observe that: 

VP-iPn) 

fi + h 

It provides a L°° bounds on ^(p„)-^^=^ since f3i + ttjt — -r^ > 0. 
The third term gives: 

g{Pn)p{Pn)DUn = 7p{pn)pt y/ pipn)^^ PiPn)DUn, (5.91) 

By p. op we prove that ip{pn)Pn y/p{pn) is bounded in L^(L°°(M^) and it yields that 
ip{Pn)Pn yj p{Pn)\/ p{Pn)Dun is Uniformly bounded in Lj,(L^(M^)). Next we have: 

'^'^g{Pn) ■ p{Pn)DUn = 2tp {pn)p^^Vpn ■ p{pn)DUn 

^,,, Jn~'^^/Kp^) p{pn)^ r^-^j. (^-^^^ 

+ 2i){pn) ttH l^^Pn ■ ^/p{Pn)DUn. 

P [Pn) yJPn 

The first term is easily bounded in Ly(L"'^(M )) since the support of V' is included in 
C(0, 1,2). In order to deal with the second term we observe that via (jl.Sp . (|1.6p and 
dLlD: 



/3i- 



WKPn) rh-. 1 = IV'(Pn) ■ ' 






It implies that V(/o«) ^"' J.^ y^ is bounded in L5?(L~(R^)) since /^i-Hw-lAf^O 

and we deduce that the second term is bounded in L|,(L^(]R^)). The two last term in 
(|5.74p are similar to treat. 

dSSSD, dSjO]), (lOB and ([5:92|) implies that ^tt;^ is uniformly bounded in L^W-^'^iK)) 
for any compact K. D 

Step 2: Convergence of ^fp^Un and Pn^^n 
In the sequel we shall define h(^p) as follows: 

M/^) = V'(p)^ + (l-V'(p))^/p if 1^1 > 2, 

\/P (5.93) 

h{p) = Tp{p)p^^+^ + (1 - V(/o))/+^ if < z^i < 2. 
Here ip, (3 and /3i verify the hypothesis of lemma [5j 
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Lemma 2 The quantity ^Jp^Un strongly converges in Lf^^{{0,T) x Q,) to j^ (which is 
null when v = 0). 

In particular, we have v{t,x) = a.e on {p{t,x) = 0} and there exists a function u{t,x) 
such that v{t,x) = y/p{t,x)h{p){t,x)u{t,x) and: 

y/p^Un -^ y/pu strongly in Lf^^{{0,T) x M^), 
PnUn -^ pu strongly in L^(0,T; L;^p^(M^)). 



Remark 26 Here u is not uniquely defined on the vacuum set {p{t,x) = 0}. We will 
set u = on the vacuum set {p{t,x) = 0}. 

Proof: Since ^J'" -. is uniformly bounded in L°°(0,T;L^(M )), Fatou's lemma implies 

that: 

v^ 
liminf ".„ (ix < +oo. 

HPnY 

We deduce that v{t, x) = a.e. in {p{t, x) = 0} since h{p) = when p = 0. We can define 
the limit velocity by u{t, x) with u{t, x) = r-,^ ^^\{^tl,u ^\\ when p{t, x) 7^ and u{t, x) = 
on {p{t, x) =0}. In particular this last point implies that there is no concentration effect 
of PnUn ® Un on the set {p = 0}. And for all t > 0: 

dx = p{t,x)\u{t,x)\ dx < +(X). 



Irn h'^ip)it,x) Jf,N 

Furthermore applying the Fatou's lemma once more, we obtain: 

/ p\u\ ln(l + |n| )dx < I liminf /9„|u„| ln(l + |u„| )dx 

< liminf / Pn\un\ ln(l + |un| )dx, 

which yields p\u\'^\n{l + jup) E L°^(0,r; L1(M^)). 

Next, since Vn and pn converge almost everywhere, we know that in {p[t,x) 7^ 0}, 

■s/]hiun = fc/"" N converges almost everywhere to ^u = jj-^. In particular it implies 

that: 

^/p^iUnl{{\u^\<M}n{p>Q} ^ \fpul{\u\<M} almost everywhere. 

(5.94) 
y/p^Unl{{\u^\<M}n{p=Q} < M ^/p^ -> almost everywhere. 

Following the argument of the proof of the lemma [1] for any compact K we have: 



\\/K.Un - y/pu\ dxdt < I \\/P^Unl{\u„\<M} - ^/pul{\u\<M}\ dxdt 

{0,T)xK J{0,T)xK 

+ 2 \./p^Unl{\u„\>M}\'^dxdt + 2 \^ul{iu\>M}\'^dxdt, 

J{0,T)xK J{0,T)xK 

(5.95) 
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Let us deal with the first term of ()5.95p we have then: 

l\/P^^nl{|«„|<A/} - Vpul{iui<M}\'^dxdt < 

{0,T)xK 

/ \VP^Unl{\u„\<M} - Vpu'i^{\u\<M}\'^dxdt (5.96) 

J{{0,T)xK}n{,/]^<c} 

J{{0,T)xK}n{y^>c} 



The first term of ()5.96p converges to when n goes to +00 via the theorem of dominated 

convergence. 

Now let us recall that via the inequality ()1.9p we have: 

Pn ^ ™ < — —— when pn > 1. 
VP^ 

Since ^^plJ is uniformly bounded in L^(L^(R^)) for A^ = 3 (see the lemma[3]) we deduce 

that ^/fhi is uniformly bounded in L^ [L'^'^'^'^ {M.^ )) . It allows to deal with the second 

term of ()5.96p for A^ = 3 by Holder's inequality and Tchebytchev lemma with c going to 

infinity. 

For N = 2 using the inequality (|1.9p and the fact that ^^^' is uniformly bounded in 



L^(L|'q^(R )) for any p > 1 when N = 2 (see the lemma[3]), we deal similarly with the 
second term of ()5.95p via Holder's inequality and Tchebytchev lemma. 
Finally, we take advantage of the gain of integrability on the velocity furnished by the 
entropy (|3.38p : 

/ \^/p^Unl{\u„\>M}\'^dxdt < / p„|u„|^ln(l + \un\'^)dxdt. 

Similarly we have: 

/ \VPul{iu\>M}\'^dxdt < i^^Q , j^,p^ / Pl^i|^ln(l + \u\'^)dxdt. 

Combining all the previous estimate, it yields: 

/C 
\^/p^Un - ^/pu\'^dxdt < 
ln(l + M^) 



for all M > 0, and the lemma follows by taking M — )• +00. 

• Let us prove now the strong convergence of pnUn to pu. Since y/p^Un converges 
strongly in Lf^^{[0,T) x M^) to ^/pu and that via the lemma [1] y^p^ converges also 
in Lf^^{{0,T) X R^) to ^ it implies that p„u„ converges strongly in Lj^^{{0,T) x M^) 
to pu. D 
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Step 3: Convergence of the diffusion terms 

Lemma 3 We have the convergence in distribution sense up to subsequence for any 
T>0: 

Kpn)Vun -^ Kp)Vu m V'{{0,T) X M^), 

KpnY^Un -^ KpYVu in V'{{0, T) X M^), 

and: 

\{pn)d\-vun -^ \{p)d\-vu in V'{{0,T) x M^). 

Proof: Let (/> be a test function, then: 

p{Pn)^Un(j)dxdt = — p{pn)Un'^ (pdxdt + / Un"^ p{Pn)<P dxdt 

= - — -^ y/p^UnV 4) dxdt + / JJh^Un — r^V Pn(t> dxdt. (5.97) 

J VP^ J y/Pn 

= - " \fpw^-nS4> dxdt + - / ^/p^UnVf{pn)4>dxdt. 

J ^JPn 2 J 

From lemma[3]in the appendix we know that ^-^^=d- is uniformly bounded in L°° (0, T; L^^^(]R )). 



Furthermore via the inequahty (|1.9p and the convergence almost everywhere from /)„ to 
p we know that ti^^SsX converges almost everywhere to ^-^ (defined to be zero on the 



vacuum set). Therefore by the lemma [TJ it converges strongly in L^^^((0,T) x M^) to 



f^- This point is enough to prove the convergence of the first term as ^fp^Un converges 
strongly. 

Next since Vf{pn) is bounded in L°°(0,T;L^(M^)), up to a subsequence Vf{pn) con- 
verges weakly to v in Lf^^{{0,T) x M ). In addition by Sobolev embedding we know that 
f{pn) is bounded in L°°(0,T; L^^^(R^))). Since f{pn) converges almost everywhere (/ is 
a continuous function) to f{p), it converges strongly in Lf^^{{0,T) x M^) by using the 
lemma [Ij It follows that: 

V/(/5n) ^ Vf{p) lL((0, T) X R^) - weak. 

It concludes the proof for the second term of (|5.97p . 

A similar argument holds for p{pnfVun and A(p„)divn„ inasmuch as we have \X{p)\ < 

Cp{p) and |A (/9)| < Cp (p) via (jl.6p and the remark [51 

Global existence ■when uq = —Vip{pQ) when p{p) = pp" with a > 1 — j^ 

The first thing to observe is that p{p) = pp"' with a > 1 — -^ verifies the hypothesis of 
theorem 11.21 for the stability of the global weak solution and in particular ()1.6p .Now it 
suffices to construct a sequence of global regular solution {pn,Un) verifying the hypothesis 
of the theorem II. 2| in particular the uniform bound via the entropy ()3.36p , p.37p and 
(|3.38p and the following properties: 



Pq converges strongly to po in L^{ 

PqUq converges strongly to po^o in L^(^'^] 
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Let K a function belonging in the Schwarz space 5(M ) with k > and J^n ndx = 1. 
We define Kn by: 

Kn = n K(n-). 

Let us take for example k{x) = (27r) 2 e 2 . By using the theorem 11.11 and by setting 
pn = pQ ^ h with /o continuous and belonging in L^(M^) n L°°(M^) n W'^'^{R^) n 
W^'^{M. ) and strictly positive and po continuous we know that there exists a global 
regular solution of the system ()1.24p with initial data (pg, — V(/3(/9q)). Indeed p^ verifies 
the hypothesis of theorem 11.21 since p^ is continuous in L^(M ) and strictly positive. We 
also observe that p^ converges strongly to po in ^^(K^)- Let us deal now with a more 
simple case when we assume that po belongs also in L°°(M^) n W^'^{M.'^) D W^'^{R^). 
Next we have: 

= a{po + ^)"-i Vpo + -{po + ^)°-tv/o. ^^'^^^ 

n n n 

Let us distinguish two cases a > | and 1 — -^a < |. 

• a > |. In this case since po belongs in H^iM.^) nL°°(M^) we deduce that V/(/Oo) and 
Y^Po'^o ^^^ uniformly bounded in L°°(L^(]R^)). Indeed we have: 

\^M^iPo)^Po\ < (« + 1)(II/'o||l-{r^) + II/o|Il-(r^))""^(|Vpo| + iv/ol), 



which implies the previous statement. 

Let us prove thatY^p^(t(Q) p is uniformly bounded in L^(IR ) for a p > 1 large enough. 
It will be then sufficient in order to show that -y//j^UQy^ln(l + IuqP) is uniformly bounded 
in L^(M^). We have then for p large enough (1 -\ — )(a — 2) + 2 > and: 

(5.99) 



Since we have assumed that po ^-nd /o are in VF"'^'^(]R ) n W^''^{W ) it shows that for p 

large enough, y^/'o^C^o) ^ i^ bounded in L^(M^) which implies the desired result. 
Let us finished by proving the strong convergence of PqUq to pqUq. We have then: 



|p>S - Pouo\ = |a[(po + ^)"-' - Po""']Vpo + ^(Po + ^)°-'V/o|. (5.100) 

Since po and /o are bounded in L°°(]R^) we deduce that: 

\\{P0 H )"~ - Po IIl°°(M^) — ^n->+oo 0. 



n 



This last inequality and Holder's inequality show in particular that [(po + „ )°~^ — 
(po)"~^]Vpo converges strongly in L^(M^) to 0. The second term of (|5.100p is easy to 
treat. It gives that PqUq converges strongly to pouo- 
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In particular via the first part of the theorem 11.21 it imphes that {pn,Un) converge in 
distribution sense to a global weak solution [p, u) . 

• 1 — jj<a<^. In this case we shall assume moreover that y^|V(/7(/9o)| ^ and 
^/fo\y^{fo)\ ^ are bounded in L^(M^) for p large enough. We have via (j5.98|) : 

\Vp^^'ip^)^P^\ < apo'Vpol +ani-"/r'|V/o|. (^-lOl) 

Since we have a > ^ it implies that y/J^(p {Po)^Po is uniformly bounded in L^(M^) by 
using the fact that ^/po'V(p{po) and \/7oVc/5(/o) are bounded in L^(R^). It implies that 
V/(po) and y^u^ are uniformly bounded in L°°(L^(M^)). Next by (j5.99p we have: 

< c{pSf^'^^^-'^^H\Vpo\'^'^ + -^WM'^h 

n^p (5.102) 

<C{pq \Vpo\ p+n"- p"- ' ^ pfo |V/o| p, 

< C(po \Vpo\ p+n^2 p^ 7o |V/o| p. 

Since q > 2 by choosing p large enough we obtain that n 2 '^ p^ '^' which is uniformly 

bounded. By the fact that Y^|Vc/9(po)| ^ and \/7o|V99(/o)| ^ are bounded in L^(M^), 

it implies that \f(^{u^) ^^ is bounded in L^(IR^) and so that pQ|nQpln(l + lugP) is 
uniformly bounded in L^(M^). 

Finally by (jS.lOOp we are going to prove that /OqUq converges strongly to pi^UQ. We start 
by remarking that when a < 1: 

i[(po + ^r-^ - pr']vpoi < 2/.o"-i|vpoi. 

We deduce by the theorem of dominated convergence that [(po + „ )°~^ ~ P%~ ^ P^ 
converges strongly to in L^(R ). The second term on the right hand side of ()5.98p goes 
also trivially to 0. It the case where 1 < a < | it suffices to observe that: 

ll|[(PoH )"~ -Po~ ]|Il°°(rjv) -^n^+00 0. 

It implies in particular since V/jq belongs in VF^'-'^(M^) that [(po + „ )"^"'^ — Po~^]Vpo 
converges strongly to in L^(M ). It achieves the proof of the strong convergence of 
PqUq to po^o- Finally it implies that {pn,Un) converge also in the case 1 — -^<a<|to 
a global weak solution (p, u) of the system (|1.24p . 

We have previously proved the existence of a global weak solution (p, u) of the system 
p.24p by assuming extra conditions on the initial density po, typically po belonging in 

L~(M^)n-H'HM^)nT^^'i(M^), po continuous and ^\Vip{po)\^^p bounded in L'^{R^) 

for p large enough. Let us deal with the general case, it suffices by a second approxi- 

i+- 
mation on the initial data (pq, ^/Pn{uQ) ^) with p large enough to pass to the limit by 
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using the first part of the theorem 11.21 More precisely by a convolution approximation 
we choose p^ = po*Kn which belongs in L°°(M^)ni?^(M-^)nM^-^'^(M'^) and is continuous 

i+i i+- 

and we set ^/fh^{uQ) p = {y/pu^ ^) * k„. 

Let us now describe the form of the solution (/?, u) , in particular we are interested in 
proving that p is also the unique global strong solution of the system I1.12p (for a such 
result, we refer to the theorem I2.4p . To do this we recall that our first approximation 
(pmUn) is solution of (J1.12p with a initial data strictly positive in L^(M^) and contin- 
uous. Let us recall that from theorem 12.41 the porous media equation verifies a crucial 
property which ensures the uniqueness the L^ contraction principle. Let us apply this 
property to the sequence (/9n)neN, we have then for all n, ?tt, G N : 

\\Pn{t, •) - Pm{t, ■)\\li^(rn^ < Wpq - p'^WL^RNy 

Since we know that Pq converge to po in i^(M^) it implies that (/5n)neN is a Cauchy se- 
quence in C([0,T],L^(]R^)) for any T > which implies that (pn)neN converges strongly 
to pi in C([0,T],L^(]R )). But since via the first part of the theorem ll.2l we know that p„ 
converges also strongly in C([0,T], L[g^(R^)) to p, it implies that p = pi. Furthermore by 
the definition of the L^ solution of the porous media equation (indeed these last one are 
defined as limit of energy global weak solution after a regularization of the initial data po 
where we use the fundamental L^ contraction principle which ensures the uniqueness of 
a such process, we called this the limit solutions in the L^ setting. We refer to the proof 
of the theorem 12.41 where we explained precisely how are defined the L^ solutions. For 
more details on the theory the reader can also consult the chapter 6 and 9 of the excellent 
book of Vazquez [33]), we know that pi is the unique solution of porous media equation 
with Pq = Pq. It proves that p = pi \s the global unique solution of ()1.12p with initial 
data Pq which belongs in L^(M^). We proceed similarly for the second approximation 
{pn,Un) since each time the approximated sequel verify the L^ contraction principle of 
the porous media equation. It concludes the proof of the theorem 11.21 D 

5.1 Proof of Theorem 11.31 and corollary [1] 

5.2 Proof of theorem 11.31 

We are now going to prove that if we have some global weak solution {p^,Ue) for the 
system (|1.23p in the sense of the definition [?] (or see |29J). then these global weak 
solution converge in distribution sense to a quasi-solution (/?, u) with initial data {po,uo). 
To prove this, it will suffice to use the same compactness argument than in the previous 
section; except that we shall deal with the pressure term and that the entropy ()3.43p 
is quite more complicated since there is a reminder term to deal with. Via the entropy 
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(fXiTD . Srm we have: 

\pJuJt,x)\'^(t,x) -\ -p2]dx+ / / fi(p^)\DuJ'^dxdt 

7 - 1 ./n JwN 



e '■' 



7 - i Jo JRN 

+ I / X{p^)\dwu^\'^dxdt < / [poho|^(2;) H 7Po]^^- 

oJm h^ 7-1 ^^^^^3^ 



[/9e|n,(t,j;)|2 + p,|Vv9(/9e)P(t,x)]dx + e / / Vipip,)-Vp2dxdt 

Jo JR^ 

< C7( / (po|noP(x) + />o|Vvp(/)o)P(x) + ^--/)2(x)) dx). 

jRr^ 7 - 1 ' 

Lemma 6 We are going to distinguish two cases. 

2. 

• 1^1 > 2 y/epe is uniformly bounded in L^(L®(M )) for any T > when N = 3 and 
in Ll{Li{R^)) for any T > and q > 2 when N = 2. 

• < z^i < 2 

1. epe ^ ^^ is uniformly bounded in Lj,{L^{E.^)) for any T > when N = 3. 

2. epT^ ^ is bounded in any L^{LP{R^)) with p € [2, +00 [ when N = 2. 



Proof: The inequality ()5.103p insures a uniform bound of -v/eA/(/3'(/>e)/52 Vpe in L^((0,T), L^(IR^)) 



for any T > 0. Let us evaluate -J (f' {pe)p7 Vp^, we have using p.6p and (jl.Sp : 



,7-1t7. -./M££)±MM.7~3^.,, 1,,.. . 7-3 



^'{p,)pr'Vp, = \l _yrii_J^l!^p'y~^ + 2(1 - -)/u(/,,)p2-^ Vp„ 



>^/2(1-1 + ^Hp,)p2-' Vp, 



Via (11.91) we deduce that: 



' "'' IVpei < c\\/^'{p,)p7-'Vp,\ \/p, > 1, 



,2 -'- 2Ar + 4iV|Y7^l^/^L/,^'('^ ^^7-l 



Let set -0 a C'o° function such that ■;/' = 1 on B{0, 1) and V' = on "^-6(0, 2). It implies 
that since suppV' is included in the shell C(0, 1,2): 

• ^V(/Oe)Vp|~^"^'^ is uniformly bounded in L|,(L2(M^)) for any T > 0. 

• \/e(l - ij{p,))Vpf~^^^ is uniformly bounded in L|,(L2(M^)) for any T > 0. 

Next by (|5.103p we know that e^ p^ is uniformly bounded in L^{L'^(W )) for any T > 

which implies that y/epi is uniformly bounded in L^(L^(M^)) for any T > 0. Let us 
deal with two different cases. 
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i/i > 2 

In this case we have : ~27v ~'~ iw' — ^ ^^ imphes that \fe{\ — ip{pe))'Vpe is uniformly 

bounded in L^(L^(R )) for any T > 0. And since e^ pj is uniformly bounded in 

L!^(L^(IR^)) for any T > 0, we deduce that e^{l — tlj{pe))pl is uniformly bounded 

in Ly (L^(IR^)) for any T > 0. We deduce that e2 (1 — ij;(p^))pi is uniformly bounded in 

LJ,{H\R^)) for any T > 0. Easily since e^tp{pe)pl is uniformly bounded in L^{L\R^))n 

LJ^(L°^(IR^)), we deduce by Sobolev embedding that e^ pf is uniformly bounded in 
Lf,(L6(M^)) when TV = 3 and in Lf,(L«(M^)) for any q>2. 

0<ui<2 

Let us distinguish the case N = 2 and A^ = 3. 



•N = 3 

We know that ^/e{l - tlj{pe))Vpi~^^*^ is uniformly bounded in L'^{L'\R^)) for any 
T > 0. Furthermore since ^/epe is uniformly bounded in L^(L^(M^)) we deduce that 

ep (1 - V'(pe))p?"^^^ is uniformly bounded in L^{LP{R^)) with p(i - 2F + i^) = 2 
with p > 2 (indeed it is possible because 2 ~ Jn ~^ In ^ ^) ■ I^ideed we have: 

1 2_J_ + ilL 1 2: 

Since p > 2 it implies that ^/e{l-tp{p^))pJ~'^^^ is uniformly bounded in L^{LP{R^)) 

and also in L^(L^(IR^)) because \/e(l — ip{pe))pi '^'^ ^"^ is strictly positive only a set of 
finite measure (it is a direct consequence of the Tchebytchev lemma) . We have shown that 

^(l-V(/?e))/?f "^"^^ is uniformly bounded in L'^{H'^{R^)) for any T > 0. By Sobolev 

embedding we deduce that e(l — ip{pe))pe ^ ^^ is uniformly bounded in L^(L^(IR )) 

for any T > 0. 

Since we know that ij{pe)pe is uniformly bounded in L^(Li(M^))nL5?(L°^(M^)) we de- 
1 . t^i 

duce that etlj{pe)pe ^ ^'^ is uniformly bounded in L;^(L'^(R )) for any T > because if 

7 — ;^ + ^ > 1 this is obvious by interpolation. In the other case etlj{pe)pe ^ ^^ is uni- 
formly bounded in L5?(LP(M^)) withp(7-^ + ^) = 1. Andwehavep= __J ^ < 2, 

we conclude also by interpolation in order to prove that eip{p^)pe ^ ^^ is uniformly 
bounded in L5?(L3(mA^)). 

Finally we obtain that ep^ ^ ' ^^ jg uniformly bounded in L^{L'^(R^^)) for any T > 0. 



^T-w + Tn : ;r l_. 1 J„^l :„ rl/T3{m,N\ 



• N = 2 

Similarly we obtain that \/e(l — tp{pe))pe '^'^ "'^ is uniformly bounded in L'^{H^{R )) 

X 1 I "1 

for any T > 0. It provides a uniform bound for -v/e(l — 'ip{pe))pe '^^ "'^ in any 

L|,(L'?(M^)) with q G [2, +cx)[. It means that €{l-7p{pe))pT^^'^ is in any L\,{LP{R^)) 
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with p G [l,+oo[. Since i^{pe)pe is uniformly bounded in L^{P{R^)) D L5?(L~(R^)) 

we deduce that e'ip{pe)pe '^ ^'^ is bounded in any L;^(LP(M^)) with p G [l,+oo[ if 
7 - i + ^ > 1 and in any L^(LP(M^)) with [pi, +oo[ where ^1(7 - ^ + ^) = 1 in the 
other case where pi < 2. It imphes that: 

• e/>r^^^ is bounded in any L].{LP{m^)) with p G [1, +cx)[ if 7 - ^ + ^ > 1. 

• epT^ ^ is bounded in any L^(LP(M^)) with pG [pi,+oo[ with pi(7-^ + ^) = 
1 in other case. 

It achieves the proof of the lemma. D 

In the following lemma we are going to get uniform estimate on the pressure p'J. 
Lemma 4 Let us distinguish two cases: 

• i^i > 2. The pressure ep2 is bounded mLs ((0, T) xM ) when N = 3 and L^{{0,T) x 
M^) for all r G [1, 2[ when N = 2. 

• < ui < 2. The pressure ep] is bounded in L^{L'^^{'R )) with r\ = 2 — a/^_^^ \ 
when N = 3. 

Proof: 

• z/i > 2: 

7 

We have seen in the lemma [6] that when N = 2, ^p^ is bounded in L^(0,r;L^(M )) 
for any q>2. We deduce that epj is bounded in Li(0,r; LP(M^)) n L°°(Li(M^)) for 
all p G [1, +00 [, hence by interpolation epl is bounded in -L'~((0, T) x M ) for all r G [1, 2[. 

When A^ = 3, by Sobolev embedding we only obtain that \/epe bounded in L^(0, T; L^(R^)) 
which gives that epl is uniformly bounded in L^(0,r;L^(M^)). By Holder inequality we 
have: 



2 3 

5 1 1 T 1 1 5 



I'YII ..^ll'YII'i ll'YII'^ 

I'^Ke II r §//^rn^.,,™^f — ll'^Pe llr,oo/-nT-f lCIRJV^^II'^/^^ " ' 



Hence epj is bounded in L3((0,r) x M^). 
• < z^i < 2: 

When A^ = 3 we know via the lemma [6] that epe ^ ^^ is bounded in L^(L^(M )). 
We have in particular that: 

epl = {epi ^^^^)p^ ^^ . 

Via lemma[3]we have seen that p^ is uniformly bounded in L^(L-'^+''i(M^)) when A^ = 3. 
We define p such that: 

,1 i^i N /I 1^1 s 6(1 + z^i) 
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By Holder's inequality we obtain with — = 3 "I — that: 



\^Pl\\L\{Ln 



\ep. 






1 "^l 

\\ „N 2N II 

Li,(L3(KiV))||P£ \\l^(Lp 



Now since epl is bounded in L\,{L^'^ {R.^ )) n L^(L^(]R^)) we have by interpolation that 
epl is bounded in Uj,{L'i {M.^ )) with: 



r 
1 _ 

K q pi 



+ 1 



It implies the following relation ^ + ^(1 — ^) = 1. In particular we obtain that epj is 
bounded in U^^ {U-^ {R^ )) with n = 2 - ^^ 



3(1+1^1) • 



D 



We are going finally to prove that y^|ne|(ln(l+|tie| )2 is uniformly bounded in L^(L (M )) 



Lemma 5 y^|ne|(ln(l + juel )2 is uniformly hounded in L!^(L (R )) for the following 



situations: 

• 1^1 > 2 and: 



< z^i < 2 and: 



5 V'} ui 



5 ^^^^ (4-..)(l+^0 .^^^3_ 
D 12 2 — z^i 

6 12 ' 6 12 ■' 

1 + |<7 ^/iV = 2. 



Proof: Let us come back to the inequality (|3.43p . we have \/6 S (0, 2): 

1 + k.P '•^ 



/ p, i^ln(l + |neP)(t,x)dx + i/ / / /i(p,)(l+ln(l + 

JR^ ^ JO JR^ 



n,r))|L'iier(t,x)(ixdt 



<C/ / fi{p,)\Vu,\^{t,x)dxdt + Cse^ (/ (^^V^)^*^^)^*- 



27-f 



(5.104) 
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Let us deal with this inequality in order to prove the uniform bound of /O^ — 'j^ ln(l + 
In^p) in L^(L^(IR )) for any T > 0. To do this we shall estimate the right hand side of 
(|5.1U4p . Using the inequality (|1.9p we obtain that: 

^1__ < _p27-2-l+jv-2iv Vp < 1. 

nip) C 

Let us distinguish two cases when z^i > 2 and when < z^i < 2. 

• z/i > 2 and TV = 3 

From the lemma HI we know that epl is uniformly bounded in L3((0,T) x M^) Pi 

L^(L^(M^)) for A^ = 3 which implies that it exists C > such that: 

et / / p~fdxdt<C. (5.106) 

JO iR'v 

27- 1 ^_ 

In particular it implies that for 5 small enough e^{^{p^>i\—f — r)^~* is uniformly bounded 
in L^((0,T) X M ) under the conditions that: 

Indeed by Tchebytchev lemma we have: 

\{X, \pe{t,x)\ > 1}\ < IIpoIIli - 

We choose p such that p{2'j — 2 ~ ^ + ^ ~ Jn^q^ ~ 1^ with 6 small enough. By Holder's 
inequality we have: 



6^ 



/ / l{p.>i}P^'^"^-iH-|^)2^da;dt<e2 /" (/ ptdx)-^\\po\\l,.^^.dt, 

<e^-^piT\\po\\mm))^i[ I elphdxdt)l. 

Jo iR^ 






1 _ 1^1 -l 2 



2 'AT 2N' 2-S 



IS 



with ^ = 1 — i and p > 1. It implies by ()5.106p that e^l|p^>i}p 
uniformly bounded in L^((0,T) x M^) under the hypothesis (|5.1U7p . 

9 (2~; ^ 1 I "*" "^1 •) 2 

Let us deal now with the term e l|p^<i}p^ ^ ^ '^ 2n>2-s ^ \x suffices to assume that: 

27-l + --^>l ^ 7>- + -. (5.108) 

^ N 2N ' 6 12 ^ ' 



Indeed we know that l{p,<i}Pe is bounded in Llf (^^(M^)) nL^(L°°(M^)) which insures 
the uniform bound of l{p,<i}P^^^"2-i+iv-2iv)2^ in L5?(L^(]R^)) by interpolation. This 
achieves the proof of the case A^ = 3. 
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The last situation consists when A^ = 3 in using the lemmaOwhen 27—1+-^ — ^ < 1+-^ 
which is equivalent to 7 < | + ^1^1. Indeed via the lemma [3] we know that pe ^^ is 
uniformly bounded in L^(L^(M^)). It implies that e^l{p^>i}P^^^"2-i+iv-2iv)23i jg ,^^1- 
formly bounded in L^(L^(R^)) when 27-l + i-^<l + ^. In order to bounde 
e l|p^<i}P^ '^2 AT 2N>2-s -^ve use the same hypothesis than in the previous case. 

• i/i > 2, TV = 2 

In this case, the situation is quite simple, indeed we know via the lemma [3] that p^ is 
bounded in L'^{L'i(S.^)) for any q > 1. In particular it implies that ^{p^>i}P ^~^ '"iv~2iv 
is bounded in L^(L^(IR^)) without any specific condition. However we shall require a hy- 
pothesis for dealing with the term lip^<iip '^~^ ^iv~2iv which is similar to the previous 



• 0<i/i<2, iV = 2 

The proof in this situation is exactly the same than in the previous case by using the 
lemma El We need only: 

27-l + ^-^>l ^ 7>4 + f. (5.110) 

• 0<i/i<2, iV = 3 

Via the lemma[3]we have seen that ep'J is bounded in L^ {U"^ (M-^)) with ri = 2 — g/]^'^^ \ ■ 
By using exactly the same argument than in the previous case we have two possibility to 
bound the last term on the right hand side of ()5.104p : 



or 



27- 


2 Z.1 
~ 3 6 


<'-.(iZ,)^ 


27- 


-4 


5 1^2 5 Ivi 

_ -1 < 7- H 

6 12 '6 12 



)7 ^ 7<2 + ^. 



(5.111) 



(5.112) 
It achieves the proof of the lemma [5j D 

We have now proved that p^ — ^^ln(l + lu^p) is uniformly bounded in Lf^J^L^{M.^)). 
We can then pass to the limit when e goes to 0, more precisely by using lemmas [TJ [2] 
we show that pei PeUe and ^p^pU^ yj p^-pU^ converges in distribution sense to p, pu and 
^fpu ® y^M and lemma E] give us the convergence in distribution sense of the diffusion 
term. Furthermore the lemmas [1] and [2] give us the following desired strong convergence: 

• pe converges strongly to p in C([0,T], L^^"(IR )) for a small enough when A^ = 3. 

• p^ converges strongly to p in C([0,T],L^^^(M^)) for any p > 1 when N = 2. 

• y/plUe converges strongly to ^fpu in L^^^((0, T) x R^)) for any T > . 
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It remains only to deal with the term e\J pj and to prove that it converges in distribution 
sense to 0. 

Lemma 6 Let us distinguish two cases: 

• When N = 3 

e°'p2 converges strongly to in L^ {Lj^^CK^ )) for any a > when e goes to for: 

7 < 1 + i/i. 

e"p2 converges strongly to in L3~'^((0,r) x M )) for i^i > 2 and for any a > 
small enough when e goes to 0. 

e"pl converges strongly to in Lj^"(L'"i~"(R )) for < z^i < 2 and for any q > 
small enough when e goes to 0. 

• When N = 2 

e°^pl converges strongly to in L"^ {L\^^{EJ^ )) for any a > when e goes to 0. 

Proof: When iV = 2 we know via the lemma [3] that p^ is bounded in L^(L'^(]R^)) for 
any q > 1. It implies trivially that e"p7 converges strongly to in L^{LJ^^{W )) for any 
a > when e goes to 0. 

When A^ = 3 we are in a similar situation when 7 < l + i^i via the lemma [3l If i^i > 2 we 
have seen in the lemma [5] that ep] is uniformly bounded in L3((0,T) x M^)), combining 
this result with the fact that ep] is uniformly bounded in L^(L^(M^)) and an interpo- 
lation argument we obtain the result that we wish. In the case < I'l < 2 we apply 
a similar argument with ri by using the lemma HI It concludes the proof of the lemma. D 

Furthermore in the same spirit by using exactly the same arguments than in the proof 
of theorem 11.21 and the previous estimate on the pressure we prove also the stability of 
the global weak solutions for the system (|1.1|) . 

5.3 Proof of the corollary [1] 

From the previous theorem we know that (pe, u^) converges to a quasi solution {p, u) and 
that />e converges strongly to p in C([0,T],Lj^^^(M^)). In particular when p{p) = pp'^ 
if we assume that there exists a unique global quasi solution we know via the theorem 
11.21 that this quasi solution verifies the porous media or the fast diffusion equation in 
function of a inasmuch as p is solution of (J1.12p . 

In particular when we assume that the initial density po has a compact support, it implies 
that when a > 1 the support of the density p remains bounded along the time. Indeed 
it consists merely of not ing that we can find a delayed Barrenblatt solution centered for 
instance at that lies on the top of po, it means: 



0<poix)<Um{T,x) VxG 



t,N 



with m large enough and r > 0. By theorem 12.41 and the inaximum principle we know 
that: 



X . , , „^ 



0<pit,x)<it + T)-^^F{j^^-^) VxG 



57 



with F{x) = {C — 2a l''"l^)+ ^ • ^^ particular it implies an information on the expan- 
sion of the support of the solution since we observe that the support of p{t, •) is included 
in a set E{t) = CB{0, M{t + r) 2 with C > 0, M > independent of t. 
Let us prove now that p^ converges strongly to p in C([0,T],L^(R^)). We know for the 
moment that p^ converges strongly to p in C([0,T],L^^^^(]R^)), it suffices to consider K 
a compact set large enough such that for any t G [0, T] we have supp/9(t, •) C K, we have 
then: 

||p(t,-)-Pe(t,-)llLl{/r)^e^0 0. (5.113) 

Now we have by conservation of the mass and the fact that pit, •) = in K'^ for any 

te[o,r]: 

||p(t,-) - Pe{t,-)\\Li{K-) = \\Pe{t,-)\\L^{K-), /r ...n 

(5.114) 

= IIPollLi(ffi^) ~ \\P(-{ii-)\\L^{K)- 

In particular since pe converges strongly to p in C([0,T],L^(i^)), ()5.114p implies that 
||/9e(t, •)||/,i(/^) converges uniformly on [0,T] to ||/o(i, •)||li(a') when e goes to 0. But since 
IIp(^) ■)llLi(_ft') = IIPollLiriRJv-) for any t G [0,r] (indeed the support of p{t,) is completely 
included in K), it induces that \\p{t, ■) — Pe{t, •)\\l^(k<') converges uniformly on [0, T] to 
when e goes to 0. 
Finally we have by using (|5.114p : 

\\p{t, •) - p,{t, •)||l1(rJV) < \\p{t, •) - p,{t, •)||l1{X) + IIPollLi(RiV) - \\p,{t, ■)\\l\K)- 

It implies that \\p{t, •) — pe{t, OIIlitr^) converges uniformly on [0, T] to when e goes to 
and we have shown that p^ converges strongly to p in C([0, T], L-'^(M^)). In particu- 
lar it implies that for e small enough p^ is the sum of a solution with compact support 
on [0,T] and of a term of small L^ norm. In this sense we can claim that the propa- 
gation speed of the free boundary of p^ is not so far to be finite at a small L^ perturbation. 
This implies in particular by interpolation that p^ converges strongly to p in C([0, T], L^ 
for any T >, any p < 1 + a with a small enough if A^ = 3 and with p > 1 if A^ = 2. 
In particular since via the theorem 12.61 we have: 



LP 



<cr'^-||pol 



with (Tp = (jvfali)+2)p ^^'^ ^v ~ (N(a-i)+2)v ' ^^ shows that up a remainder term of 
small norm in L^, the L^ norm of p^, decrease in time for small e. It means that in some 
sense the density is subjected to a damping effect in time for the L^ norm which is very 
surprising since this effect seems purely non linear. 

Let us deal now with the time asymptotic behavior of p^. We expect that Pe{t,-) goes 
asymptotically in time to the Barrenblatt solution Um of (|1.12p of mass ||/5o||li(rjv) = w- 
We have then: 

\\Um{t, •) - Peit, OIIliCRJV) < \\Um{t, •) - p(t, OIIl^CR^) + llp(*> O " Pe{t, OlliiCR^)- 

(5.115) 
Via the theorem 12.71 we know that \\Umit, •) — pit, OlliirR^) converges asymptotically to 
when t goes to +c«. The second term converges also to when e goes to since we have 
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shown that p^ converges strongly to p in C([0,T],L^(IR )) for any T > 0. In particular 
it implies that for any a > it exists T > such that: 

||t^m(i, •) -p(i, OIIlHr^^) < a Vt>T. 
Furthermore it exists eo > such that for all e < eo we have: 

\\p{t,-)-p,{t,-)\\m^N.<a VtG[0,nr], withnGN. 



It implies that for all a > it exits T > such that for all n € N it exits eg > such 
that for all < e < eo we have: 

\\Um{t,-) -Pe(i, •)||Li(RiV) < 2a Vt G [T,nT]. 

In this sense we observe that for e small enough the solution /j^ tends to converge asymp- 
totically to a Barrenblatt solution of mass ||po||li(R^)- I— ' 
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